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I . INTRODUCTION 


This  report  describes  plotting  routines  which  were  written  to  han- 
dle data  generated  by  two  finite  element  computer  programs,  Epic-2*  and 
Epic-32,  and  create  cross-sectional  or  isometric  deformation  plots  of 
the  impact  situation  as  predicted  by  the  programs  on  a Calcomp  Plotter 
or  Tektronix  Display  Terminal.  Geometry,  in  the  Epic  programs,  is 
modeled  as  a series  of  triangular  elements  (Epic-2)  and  tetrahedral 
elements  (Epic-3)  defined  by  three  and  four  nodes  respectively.  The 
displacements  which  may  occur  as  the  impact  progresses  are  reflected 
in  the  changes  in  the  nodal  coordinates.  (See  Figures  1 and  2). 

The  nodes  which  define  each  element  and  their  present  coordinates 
constitute  the  plot  files  created  at  desired  intervals  during  the 
running  of  the  Epic  Programs.  Although  the  plotting  programs  described 
herein  were  written  to  handle  Epic  plot  files  specifically,  they  could 
readily  be  adapted  to  interface  with  other  programs  which  generate 
the  same  geometrical  data. 

These  plotting  programs  were  written  in  standard  Fortran  for  the 
Univac  1108  Computer  but  could  easily  be  utilized  on  most  other  com- 
puters of  equal  size.  It  is  assumed  that  the  potential  user  of  these 
plotting  programs  has  access  to  either  (A)  a Tektronix  Display 
Terminal  which  communicates  with  the  computer  on  which  the  Epic  program 
is  run  and  the  Tektronix  Plot  - 10  software  is  resident  (Note  changes 
to  this  in  Appendix  A for  this  program.)  or,  (B)  a Calcomp  Plotter 
wherein  the  standard  Calcomp  software  is  resident  on  the  computer  and 
automatically  creates  a file  to  drive  the  plotter. 

The  original  cross-sectional  plotting  program  was  supplied  by 
Honeywell  with  the  Epic-3  code.  It  was  modified  by  the  author  as 
discussed  in  Section  II  and  is  included  in  this  report  primarily  for 
completeness . 


II.  THE  CROSS-SECTIONAL  PLOTTING  PROGRAM 

A.  Options. 

This  plotting  program  handles  plot  files  created  by  Epic-2  or 
Epic-3  and  generates  a cross-sectional  deformation  plot  of  an  impact 
situation.  All  of  the  nodal  coordinate  connections  are  plotted  for 
Epic- 2 data  whereas  only  those  nodal  coordinate  connections  between 
nodes  for  which  Y = 0 are  plotted  for  Epic-3  data.  The  original 
plotting  program,  provided  by  Honeywell,  generated  a plot  tape  for  the 
Calcomp  plotter.  This  version  of  the  program  (listed  in  Appendix  B) 

1.  Johnson,  Gordon  R. , "Epic-2  A Computer  Program  For  Elastic-Plastic 
Impact  Calculations  in  2 Dimensions,"  Honeywell,  Inc.,  Dec.  1977. 

2.  Johnson,  Gordon  R.,  "Epic-3  A Computer  Program  For  Elastic-Plastic 
Impact  Calculations  in  3 Dimensions,"  BRL  CR  343,  July  1977. 

(AD  #A043281) 
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Figure  1.  Geometry  in  the  Epic-2  Program  (Courtesy  of 

3 

Gordon  R.  Johnson  ) 


3.  "Epic -2  A Computer  Program  for  Elastic-Plastic  Impact  Calculations 
in  2 Dimensions",  Interim  Technical  Report,  March  1977. 
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Figure  2.  Geometry  in  the  Epic -3  Program. 

(Courtesy  of  Gordon  R.  Johnson2) 


has  been  modified  to  provide  the  following  options: 

1.  handle  output  from  the  Epic-2  code  also.  (The  GPLOT 
subroutine  in  Epic-2  was  accordingly  modified  to  interface  with  this 
program.  See  Appendix  C) ; 

2.  produce  plots  on  a Tektronix  4014  Display*  Terminal,  (see 
Figure  3); 

3.  create  reflection  about  the  z axis  for  Epic-2  generated 
data,  (see  Figure  4). 

B.  Data  Format. 

Besides  the  plot  file  which  contains  the  output  generated  by  the 
Epic  programs,  the  following  input  data  is  required,  either  as  part  of 
the  runstream  or  as  a separate  file,  in  order  to  produce  cross-sectional 
deformation  plots: 

Card  image  1:  Format  (215  , 4F10.0,  611,  4X,  3A6,  A2) 

NCASE  - Case  number  (09999  for  end  of  data) 

NCYCLE  - Cycle  number 

ZMAX  - Maximum  z coordinate  included  in  plot. 

ZMIN  - Minimum  z coordinate  included  in  plot. 

XMAX  - Maximum  x coordinate  included  in  plot. 

XMIN  - Minimum  x coordinate  included  in  plot. 

IE  - '1'  if  all  elements  in  the  elastic  range  are  to  have  ’ E 1 
written  in  their  centers. 

IP  - '1'  if  all  elements  in  the  plastic  range  are  to  have  'P' 
written  in  their  centers. 

IF  - '1'  if  all  elements  which  are  failed  in  shear  and  tension 
are  to  have  ' F'  written  in  their  centers. 

ID  - '2'  if  Epic-2  produced  the  plot  file. 

IR  - '1'  if  reflection  about  z-axis  is  to  be  plotted  also 

(Epic-2  only).  (See  Figure  4.)  Note:  XMIN  will  probably 
need  to  be  changed  to  accommodate  the  reflection. 


ICAL  - '1'  if  Calcomp  tape  to  be  generated.  Otherwise  plot  will 
be  output  to  Tek  4014  screen  immediately. 


0) 

XL 

4-> 


NOTE: 


r 


i 


1.  If  the  plots  are  to  be  generated  on  a Tektronix  Display  Terminal, 
each  cycle  specified  correctly  will  be  plotted  and  a copy  automatically 
made  on  the  hard  copy  unit  (assuming  one  is  connected  to  the  terminal) . 
Press  RETURN,  the  screen  will  be  erased  and  the  next  cycle  plotted,  etc. 

2.  You  may  specify  as  many  cycles  as  desired  as  long  as  they  are 

in  ascending  order  by  cycle  number  and  the  cycles  exist  on  the  same  file. 

3.  If  you  do  not  get  anything  plotted,  check  your  output  file  (6). 
You  have  probably  specified  a case  and  cycle  number  which  are  different 
from  those  on  the  file. 

C.  Example. 

Figure  5A  is  the  runstream  for  the  Univac  1108  Computer  and  Figure 
5B,  the  input  data  which  was  used  to  generate  the  plot  shown  in  Figure  3. 


III.  THE  ISOMETRIC  PLOTTING  PROGRAM 

A.  Options. 

This  plotting  program  (listed  in  Appendix  D)  handles  output  from 
Epic-3  only  and  generates  an  isometric  deformation  plot  of  an  impact 
situation  with  the  use  of  a hidden  line  algorithm. 4 

Options  include: 

1.  creating  a reflection  of  the  nodal  coordinates  about  the 
x-axis  so  that  a complete  projectile  and  target  may  be  viewed,  assuming 
that  the  geometry  specified  is  for  half  the  impact  situation  (the  nodes 
are  restrained  in  the  y direction)  and  every  nodal  y coordinate  > 0 
(see  Figure  6 for  coordinate  system  used) ; 

2.  scaling,  translating  and  rotating*  the  scene  (specified 
in  input  parameters) ; 

3.  the  following  possible  outputs: 

a.  a plot  tape  for  future  plotting  on  the  Calcomp 
plotter  (batch  mode) , 

b.  the  scene  plotted  immediately  on  the  Tektronix  screen 
(interactive  mode) , 

c.  a plot  tape  (batch  mode)  for  future  plotting  on  the 
Tektronix  screen  (interactive  mode) . (See  Appendix  E for  the  retrieval 
runstream  and  Appendix  A for  changes  to  Tektronix  Plot-10  routines  to 
handle  this  situation) . 

T!  The  algorithm  and  basis  for  the  hidden  line  subroutine  (PLOT3D)  is 
described  in  "Algorithm  483  Masked  Three  Dimensional  Plot  Program 
With  Rotations,"  Steven  L.  Watkins,  Communications  of  the  ACM, 

Vol.  17,  Number  9,  September  1974. 

* See  Section  IIIC  for  limitation  on  rotations. 
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IASS, AX  PRPLOT.  PLOT  FILE  FROM  EPIC-3  OR  EPIC-3 
fiUSE  3 PRPLOT 

•ASG,AX  RRDUM**  CASE  AND  CVCIE  NO*  FOR  FILES  ON  FILE3 
•USE  6, RRDUM. 

•ASG,AX  R2DATA.  INPUT  DATA 
•USE  7.R2DATA. 

•PREP  RPLOT.,RLIB. 

•MAP  ,TPFf*ABS 
IN  RPL0T.2DPL0T 
IN  RPLOT.TOP 
LIB  RLIB. 

LIB  SEAPPDXTEKLIB2. 

•XQT 


NOTE:  RLIB  and  SEAPPD*TEKLIB2  contain  Tektronix  Plot-10  routines. 


Figure  5A.  Runstream  to  Generate  Figure  3. 


114  380 
09999 


*8 


-0.8 


2*8 


-2*52 


Figure  SB.  Input  Data  (in  File  R2DATA)  to  Generate  Figure  3. 
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B.  Data  Format. 


The  following  input  data  is  required,  either  as  part  of  the 
runstream  or  as  a separate  file,  in  order  to  produce  isometric  defor- 
mation plots: 

Card  image  1:  Format  (215,  4F10.0,  II,  9X,  3A6,  A2) 

NCASE  - Case  number  (09999  for  end  of  data) 

NCYCLE  - Cycle  number 

ZMAX  - Maximum  z coordinate  included  in  plot. 

ZMIN  - Minimum  z coordinate  included  in  plot. 

XMAX  - Maximum  x coordinate  included  in  plot. 

XMIN  - Minimum  x coordinate  included  in  plot. 

IOUTPT  - '1'  if  Calcomp  plot  tape  to  be  generated. 

'2'  if  plot  to  be  output  to  Tek  4014  screen  immediately. 

'3'  if  plot  tape  to  be  created  for  later  retrieval  and 
plotting  on  Tek  4014  screen. 

TITLE  - Title  to  be  written  on  plot. 

Card  image  2:  Format  (1315) 

NRING  - Number  of  rings  in  projectile  (2-5). 

NONODE  - Number  > largest  node  number  for  target. 

LNPROJ  - Highest  projectile  node  number.  (It  is  assumed  that 
the  projectile  node  numbers  are  lower  than  the  target 
node  numbers.) 

ISTARG  - Lowest  target  node  number. 

IR  - ’1'  if  reflection  about  x axis  desired.  'O'  otherwise. 

(The  reflection  can  only  be  requested  if  nodes  were 
restrained  in  the  y direction  in  Epic-3.) 

ITARGT  - '1'  if  a target  is  part  of  the  scene.  'O'  otherwise. 

NX  - the  number  of  nodes  in  the  x direction. 

NY  - the  number  of  nodes  in  the  y direction. 

NZ  - the  number  of  nodes  in  the  z direction. 
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Card  image  3:  etc: 


r 
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Scaling,  translation  and  rotation  are  specified  in  the  following 
formats.  The  order  in  which  they  are  placed  should  be  the  reverse  of 
the  order  in  which  the  transformations  are  to  be  performed,  (i.e.,  the 
first  transformation  should  be  the  last  set  of  data) . 

Scaling: 

Card  image  a:  Format  (II) 

•1' 

Card  image  b:  Format  (314) 

151  - x scaling  factor  (x  100) 

152  - y scaling  factor  (x  100) 

153  - z scaling  factor  (x  100) 

Note:  If  scaling  is  done  for  at  least  one  dimension  it  must  be 

done  for  all  three  (e.g.,  if  you  want  the  z dimension  doubled  specify 
'1',  '010001000200'). 

Translation: 

Card  image  a:  Format  (II) 

'2' 

Card  image  b:  Format  (314) 

IT1  - x translation  (x  100) 

IT2  - y translation  (x  100) 

IT3  - z translation  (x  100) 

Rotation:  Rotation  is  handled  by  rotating  one  axis  into  another 
(e.g.,  to  rotate  positively  about  y axis,  rotate  x axis  into  z axis) 

Card  image  a:  Format  (II) 

'3' 

Card  image  b:  Format  (314) 

IR1  - number  of  first  axis 
IR2  - number  of  second  axis 

ITHETA  - angle  of  rotation  (in  whole  degrees  only) . 

17 


End  of  Transformations:  If  any  transformations  are  performed 
they  must  be  followed  by  this  image  indicating  the  end  of  transformations. 

Card  image  a:  Format  (II) 

.41 


NOTE:  See  notes  1 - 3 in  Section  II. 

C.  Examples. 

Plots  generated  for  two  target  impact  situations,  one  involving 
a long  rod  penetrator,  the  other  a sphere,  and  their  runstreams  and 
input  data  are  shown  in  Figures  7-10. 

It  is  usually  more  impressive  to  see  an  isometric  plot  instead  of 
a cross-sectional  plot  of  an  impact  situation.  However,  in  cases 
where  a projectile  causes  deformation  of  a target  but  does  not  main- 
tain continuous  contact  with  the  target,  it  may  be  more  meaningful 
to  generate  a cross-sectional  plot  rather  than  an  isometric  plot. 

(Compare  Figures  3 and  11  for  which  the  same  plot  file  was  used.) 

D.  Procedure. 

The  procedure  within  the  program  is  as  follows: 

1.  Read  the  nodal  coordinates  representing  the  nodal  displacements 
of  the  projectile  and  target  at  a given  cycle  (time  step). 

2.  Create  reflection  of  these  coordinates  about  the  x-axis  if 
requested. 

3.  Perform  transformations  on  coordinates  according  to  scaling, 
translation  and  rotation  parameters  specified. 

4.  Perform  the  following  steps  for  the  projectile: 

a.  The  four  triangles,  inherently  created  when  the  four  nodal 
points  of  each  tetrahedral  element  are  joined,  are  checked  against  all 
the  other  triangles  which  make  up  the  figure.  If  duplicated,  the 
triangle  is  an  interior  triangle  and  is  removed  from  further  considera- 
tion. 

b.  Check  all  nodes  remaining  (defining  exterior  triangles) 
against  all  exterior  triangles  to  see  if  they  are  hidden  by  a triangle, 
(see  Appendix  F for  details  of  this  procedure).  If  a node  is  hidden, 
any  triangle,  of  which  it  is  a vertex,  i>s  removed  from  further  considera- 
tion. 

Almost  all  the  lines  are  now  eliminated  which  are  hidden  from  view 
for  an  individual  object,  that  is,  without  regard  to  the  relationship 
between  projectile  and  target. 
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Isometric  Plot  of  Rod  and  Target  With  Reflection  About  X Axis 


«ASG,AX  ZZPL0T2.  PLOT  FILE  FROM  EPIC-3 
AII5F  3 ZZPL0T2 

•ASG#AX  RRDUH. * CASE  AND  CYCLE  NO.  FOR  FILES  ON  FILE  3 
•USE  6/RRDUM. 

iASG,AX  RRDATA.  INPUT  DATA 
«Aii,AXRRDUnSS.  FILE  CREATED  IF  IOUTPT  • 3 


(USE  22.RDUH22. 
(PREP  RPLOT . 

(MAP  ,TPF*.ABS 
IN  RPLOT. 3DPL0T 
LIB  RPLOT. 

LIB  RLIB. 

LIB  SEAPPDITEKLIB2. 
(XQT 


(NOTE:  RLIB  and  SEAPPD*TEKLIK2  contain  Tektronix  Plot-10  Routines) 


Figure  8A.  Runstream  to  Generate  Figures  7A  and  7B. 
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Figure  8B.  Input  to  Generate  Figure  7A  (in  file  RRDATA) 
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Figure  8C.  Input  to  Generate  Figure  7B  (in  file  RRDATA) 
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Figure  10A.  Input  to  Generate  Figure  9A. 
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Figure  10B.  Input  to  Generate  Figure  9B. 
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sometric  Plot  of  Impact  Situation  Wherein  There  is  not  Continuous  Contact 


c.  Depending  on  the  angle  from  which  the  scene  is  viewed,  there 
remain  a few  lines  describing  triangles  whose  3 vertices  are  seen  but 
whose  lines  should  be  hidden  (e.g.,  connections  between  nodes  defining 
triangles  at  the  corners  of  the  target) . These  situations  are  handled 
on  an  individual  basis. 

d.  Plot  the  lines  describing  the  projectile  and,  if  there  is  a 
target  in  the  scene,  note  for  every  transformed  x value  within  the  area 
of  the  projectile,  the  highest  and  lowest  transformed  z values.  These 
values  form  the  hidden  line  mask  for  the  target,  the  area  in  which  none 
of  the  target  will  be  seen. 


5.  Perform  steps  4a-4c  for  the  target.  For  each  line  involved  in 
plotting  the  target,  compare  the  two  transformed  x values  (xA>  xfi) 

defining  it  with  the  range  cf  x values  included  in  the  mask.  There  are 
two  possibilities: 


1.  x^  and 

either  (a)  xA  > range  and  x 


vi.vuvi  j a . ' x uu^v  auu 

The  line  is  plotted  normal 


ai  c uu  toiuc  Lite  iaiigc  ui  Lite  mabh  auu 

> range  or  (b)  < range  and  xfi  < range. 


ii.  at  least  part  of  the  line  lies  within  the  range  of  the  mask. 
Starting  at  one  end  of  the  line  (x  ) check  the  corresponding  z value 
for  each  x where 


x 

x 


1 

2 


Ax 


x„=  x + Ax  and  Ax  = ± .001* 

D D - 1 

to  determine  which  segment(s)  of  the  line  lie(s)  above  or  below  the 
mask.  Then  plot  the  segment (s)  not  hidden  by  the  mask. 


E.  Limitations. 


1.  Since  the  projectile  determines  how  much  of  the  target  will  be 
seen,  the  target  can  never  hide  the  projectile  and  a positive  rotation 
from  the  z-axis  to  the  y-axis  (about  the  x-axis)  will  produce  incorrect 
results.  Similarly,  for  cases  of  significant  deformation  of  the  target, 
one  may  not  specify  reflection  about  the  x-axis  (see  Figure  13)  so  in 


*See  Figures  12A  and  12B.  There  is  a significant  difference  in  the 
straightness  of  the  lines  between  plots  where  Ax  = .01  and  Ax  = .001. 
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Section  IIIB,  Card  image  2 IR  must  = 0.  This  problem  does  not  occur 
when  there  is  no  target  in  the  scene.  This  plotting  program  should 
work  well  for  any  other  rotations  wherein,  for  6 being  the  angle  of 
rotation,  - 90°  <0<  90°. 

2.  Node  numbers  must  be  less  than  5 digits  (9999  max).  Other- 
wise the  numbering  scheme  for  the  ITRIAN  array  must  be  changed  or  a 
computer  with  a larger  word  size  must  be  used. 

3.  There  cannot  be  a total  of  more  than  2000  nodes  for  the 
projectile  and  target  together.  Otherwise  the  dimensions  of  TRANSF, 
INTEL,  NODE,  and  XYZ  arrays  must  be  changed,  (See  Section  HIE. 6). 

4.  After  the  exterior  triangles  are  eliminated,  there  cannot  be 
more  than  1600  triangles  (if  reflection  about  x axis)  or  3200  triangles 
(if  no  reflection  about  x axis)  plotted  for  either  the  projectile  or 
the  target.  Otherwise  ITRIAN  and  the  parameter  LIMIT  must  be  increased. 

5.  Since  this  program  plots  only  complete  lines  between  points 
in  view,  for  cases  of  extreme  deformation  of  the  projectile  there  are 
some  angles  at  which  a portion  of  a line  should  appear  but  will  not, 
or,  conversely,  when  a line  appears  which  should  be  hidden.  Usually 
if  the  angle  of  rotation  is  changed  slightly  (5°  - 10°)  this  situation 
will  no  longer  occur  (See  Figures  14A  and  14B) . 


6.  This  program  just  fits  a Univac  1108-65K  memory  computer 

<50  if  anv  ai-rava  arp  f»nl  arood  . mpmnrv  mnat  hfi  used. 


Figure  14B.  Isometric  Plot  of  Projectile  Shown  in  Figure  14A  Without 
Hole  (Slightly  Different  Rotation) 


IV.  SUWARY 


« 


J 


These  plotting  routines  were  developed  primarily  to  enable 
relatively  quick  and  easy  analysis  of  the  deformation  history  of  an 
impact  as  predicted  by  the  Epic  codes  by  providing  a visual  presenta- 
tion of  the  nodal  displacements.  The  key  to  the  success  of  the 
isometric  plotting  program  is,  of  course,  the  hidden  line  routine. 
Without  it,  there  would  only  be  a jumbled  pile  of  lines.  Ideally, 
one  would  employ  a more  sophisticated  hidden  line  algorithm  which 
would  handle  every  possible  situation  and  angle  of  view  in  creating 
isometric  plots  but  that  would  require  a great  deal  of  computer  time 
and  would  be  extremely  expensive.  The  techniques  used  in  this  program 
were  geared  to  the  specific  scene  arrangement  we  wanted  for  the  impact 
simulations  we  handle  and  proved  highly  successful  in  terms  of  speed 
and  cost. 

It  is  hoped  that  other  Epic  users  can  profitably  utilize  the 
programs  documented  here  to  aid  in  the  analysis  of  the  deformations 
of  their  particular  situations. 
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APPENDIX  A* 


REVISED  TEKTRONIX  PLOT  10  SUBROUTINES 
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C******** ***************************************************** ********** 

C*  PRODUCT  062-152B-02  PREVItw  ROUTINES  FOR  CALCOMP  * 

C*  RELEASE  1.1  * 

C*  C COPYRlOhT  1473  TEKTRONIX.  INC.  * 

C*  ALL  RIOhTS  RESERVtU  * 

C*  * 

C*  TEKTRONIX.  INC.  * 

C*  P.  0.  BOX  500  * 

C*  BEAVERTON.  ORtOON  4/005  * 

C*  * 

C*********************************************************************** 

c 

C BOUUNE — P 1 ItMRONlX.  INC. 


SUBhOUTINE  PIN(X.Y.N) 

REAL  A < 100) 

DATA  IFIKST.NBMAX/4B4B4J. 100/ 
IF  < IF IRST ) lu0.30U.lU0 

C*  r IRST  PASS OPEN  F ILt 

luO  IFIRST  = 0 

CALL  IFlLt (22. 'CTMP1 • > 

*f ILL  BUFFER 
2u0  WEAO(22.ENU=400) A 
Nh  = 0 

3u0  I F ( NB  .1>T.  NBMAX-3)00  TO  £00 
X = A ( No ♦ 1 ) 

Y = A(Nri*2) 

N = A(NH»J) 

NB  = NH  ♦ 3 

IF (N  .Nt.  44V) 00  TO  SUO 
AuO  N = 444 

IF  I WST  =484843 
C ENDFILE  22 

5u0  RETURN 
ENO 


Appendix  A-l.  Listing  of  Tektronix  Routine  PIN 


*Pemiaaion  was  granted  by  Tektronix,  Inc.,  May  1978  to  document  the 
following  subroutine e in  this  report. 
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PRODUCT  062-1526-02  PREVIEW  ROUTINES  FOR  LALCOMP 

RELEASE  1.1 

C COPYRIGHT  1973  TEKTRONIX*  INC. 

ALL  RIGHTS  RESERVtO 

TEKTRONIX.  INC. 

* P.  0.  BOX  500 

« BEAVERTON*  OREGON  97005 


ROUriNE— POUT TEfVTRONlX.  INC. > 

SUHROUTINt  POUT ( X * Y.N) 

REAL  A ( 1 00 ) 

DATA  IFIRST .NBMAX/9B9B93. 100/ 

IF  < IF IRST)  100.200.100 

C*F i WST  PASS OPEN  FILE 

100  IFIRST  = 0 

C CALL  OFILE (22. 'CTMPi • ) 

2o0  IF (Nb  .GE.  NBMAX-3) uO  TO  300 
A(NB*1)  = X 
A (Nb^2)  = Y 
A (NB*3)  = N 
NB  = NB  ♦ 3 

IF (N  .EG.  999 ) GO  TO  900 
BO  TO  500 
C <*UUMP  bUFFER 
300  WRITE (22) A 
Nb  = 0 
GO  TO  200 
900  «(R  I TE  ( 22 ) A 
NB  = 0 

C ENOTILE22 

b u 0 RETURN 
END 


Appendix  A-2.  Listing  of  Tektronix  Routine  POUT 
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c* 

PPOOUCT  062-1526-02  PPtVlcw  POUT1NES  FOP  LALCOMP 

c* 

PELEASt  1.1 

« 

c» 

C COPYPIGHT  ly73  TEkTKONIX.  INC. 

♦ 

c* 

ALL  pIgmTs  PEStRVtO 

« 

c* 

« 

c* 

TEKTPONlX.  INC. 

c* 

P.  0.  HUX  500 

« 

c* 

wEAVEPTON.  OPtGON  y/oos 

c* 

* 

SUwPOUTINE—PLUT TEKTPONIX.  INC. 

SUHPOuTlNt  PLOT ( X IN. Y 1N.NPEN) 

C * ThIS  POuTINE  FOLLOWS  (Ht  NOpMAL  CALCOMP  CONVENTION  OF 
C * 2=  PEN  DOWN  OP  V 1SIBLE  VECIOP 

C * 3=  PEN  UP  OP  NONVlSIttLE  VECTOP 

C * -2  = NE*  UP  IG I N VISIhLE  VECTUP 

C * -3=  NErf  OPIGIN  NONVlSIhLE  VtCTOP 

LOGlCttL  SKIP 

COMMON  /CLCUMP/  XSTaPT,YSTAPT»aFSET»YFSET»XALUM,yaCUM, 

♦ SK IP.NSKlP. IOPT.XLEN. YLEN. XF  aC. YF  AC.FAC.NhAkU 
LOGICAL  SAVE 
LOGICAL  ANUPLT 
l)AT«  SAVE/.  TPUE  ./ 

data  anopl t/. false./ 
x=x  In 
Y = Y IN 

N= I ABS ( NPEN ) 

IF  (N  .EO.  y-.b)  GO  TO  yyb 
I F ( N .Eu.  yy?) GO  TO  yy/ 

IF  ( N .EO.  yyy>  60  Tu  yyy 

IF (SKIP  .AND.  NPEN  .6T.  0)  GO  TO  POO 

NSK1P=NSKIP-1 

SKIP=NSKIP  .GT.  0 

C IF  IT  IS  UESIPEU  TO  APPLY  OFFSETS  ANO  SCALING  TO  ALL  VECTOPS 
C PC.MOVE  NtXT  STATEMENT 

IF  (N  .Nt.  12  .AND.  N .NE.  1 3 > oO  TO  UOl 
C APPLY  OFFSETS  ANO  FACTOPb 
N=N-1(J 

X=X*XF AC*XFSET 
Y=Y*YF AC* YFSET 
oni  X=X*F AC*XACUM 
Y=Y*F  aC* YALUM 

IF (SAVE)CALL  POUT (X-aACUM.Y-YaCUM.nPEN) 

IF(ANoPLT)GO  TO  100 
IF (N  .tO.  3)  GO  TO  UOJ 
IF (N  .EO.  2)  GO  TO  002 
C LOPOP  IN  CALL 
GO  TO  yoo 

C UPrtW  IF  NPtN=2.  -2.  12.  -12 

002  CALL  OPAftA(A.Y) 

GO  TO  100 

C MOVE  IF  NPEN=3.  -3.  13.  -13 

003  CALL  MOVE  A ( X » Y ) 

Appendix  a-3.  Listing  of  Tektronix  Routine  PLOT 
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C F INAL  CHtCK 

100  IF ( NPEN  . GT . 0)  GO  TO  900 
IF (SKIP)  GO  TO  900 

C tXERCISE  OPTION 

GO  T0(101*102*103tl05)  » IOPT 
C i =CmaNGE  ORIGIN  ANO  PROCEDE 

101  XACUMxX 
YACUM=Y 
GO  TO  900 

C <=CHANGE  ORIGIN  ANO  INuUIRE 
10^  XACUMsX 
YACUM=Y 
GO  TO  104 

C j=CHANGE  TO  USER  ORIGIN  ANO  INQUIRE 
103  XACUM=XSTAkT 
YACUMxYSTART 
1<M*  CONUNUt 
C **=uStR  OtFlNti)  OPTION 
C UStR  MAY  INSlRT  COUE  MtRt 
C uPTION  MOST  dE  ALTEREO  TO  PERMIT  4 


10S 

c 

99a 


900 

999 


GO  TO  900 
CHANGt  SAVE  F LAG 

SAVt=. NOT. SAVE 
GO  TO  900 
ANOPLT  = .TRUE. 

GO  TO  900 
RETURN 

if  ( save ) call  poutu.y.n) 
call  FIN1TT (0.0) 

GO  TO  900 
ENO 


! 
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* PRODUCT  062-1S26-02  PREVItrf  ROUTINES  FOK  CALCUMP 

* RELEASE  1.1 

* C COPYRIomT  1973  TE*THONiX.  INC. 

* ALL  RIGHTS  RESERVtD 
« 


TE*VTRONl X » INC. 

P.  0.  BOX  SOO 
BEAVERTON.  ORtGON  */ 00b 


RUUTINE— MAINLINE  RETRIEVAL  S YS T tM- Ttiv  TRON I X . INC. 

CALL  INITT(JO) 

CALL  PLOT (U. .0. .99b) 
luO  CALL  PLOTS ( U . 0 . 0 ) 

2u0  CALL  PIN(X.Y.N) 

IF (N  .EU.  999 ) GO  TO  9*9 
CALL  PLOT(X.Y.N) 

GO  TO  200 
C*  tND  OF  PASS 

9*9  CALL  TOuTPT ( 7 ) 

CALL  TOUTPT (7) 

CALL  hOCOPY 
CALL  TINPuT(K) 

C * A »0»  WILL  STOP  The  JUrt 
If  (K  .NE.  Hi ) GO  TO  100 
CALL  UNIT!  (0.0) 

END 


Appendix  A-4.  Listing  of  Tektronix  Mainline  Retrieval 
Routine  - RETREV 
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APPENDIX  B.  LISTING  OF  CROSS-SECTIONAL  PLOTTING  PROGRAM 


THIS  PROGRAM  PLOTS  DEFORMED  GEOMETRY  FROM  EPIC-2  OR  EPIC-3  OUTPUT. 
IT  PLOTS  THE  X-Z  PLANE  GEOMETRY  AT  Y « 0. 

THE  DATA  IS  OBTAINED  FROM  TAPE  3,  WRITTEN  IN  SUBROUTINE  GPLOT, EPIC-2 
OR  TAPE  13,  WRITTEN  IN  SUBROUTINE  GPLOT, EPIC-3 


DIMENSION 

1 NODE! 4000) , X( 4000) ,Y(4000) , Z (4000) , I F IX ( 4000) , 

* SYM( 3 ) , SYMB ( 3) , XSI ZE (4) , ZSIZEI5),  XP<6),  ZPI6J 

DIMENSION  TITLE(4),IBUF(1000) 

DATA  XX, YY/1HX,1HY/ 

DATA  BLANK/6H  / 

DATA  SYM ( 1),SYM(2),SYM(3)/  1HE,1HP,1HF/ 

DATA  LIT1,LIT2,LIT3,LIT4,LIT5/49,50,51,52,53/ 

DATA  IFIRST/O/ 

YT ITLE=YY 
XT ITLE=XX 

100  F 0RMAT(2I5,4F 10.0, 611, 4X,3A6,A2) 

101  FORMAT  (/////, 15X,23HEND  OF  TAPE  3 PLOT  DATA) 

102  FORMAT  ( 1H1 , 5X, 7HCASE  « , 1 5 , / ,5X,7HCYCLE  » , 15 , / , 5X , 7H IPLOT  «,I5) 

103  FORMAT  ( 4X, I 5 , 3E 1 5. 6 , 5X, I 5 ) 

104  FORMAT  ( 1H1 ,////, 9X , 38HE LE  NODEl  N0DE2  N0DE3  N0DE4  ICHECK,///) 

105  FORMAT  (5X.6I7) 

INITIALIZE  PLOT  AND  ESTABLISH  ORIGIN  AT  LOWER  LEFT 

REWIND  3 
REWIND  7 


READ  INPUT  DATA  FOR  PLOT 


**** 
**** 
C***« 
C***« 
C**«* 
C**** 
C**** 

c**»* 

c***« 

c**** 

c**»* 

150 


153 


IE  = l IF  ALL  ELEMENTS  IN  ELASTIC  RANGE  ARE  TO  HAVE  ' E*  WRITTEN  IN 
THEIR  CENTERS. 

IP  = 1 IF  ALL  ELEMENTS  IN  PLASTIC  RANGE  ARE  TO  HAVE  *P*  WRITTEN  IN 
THEIR  CENTERS. 

IF  * l IF  ALL  ELEMENTS  WHICH  ARE  FAILED  IN  SHEAR  AND  TENSION  ARE  TO 
HAVE  • F • WRITTEN  IN  THEIR  CENTERS. 

ID  = 2 IF  2D  VERSION  OF  EPIC  PRODUCED  PLOT  TAPE 

IR  x 1 IF  AXISYMMETRIC  PROBLEM  AND  REFLECTION  ABOUT  Z AXIS  IS  TO 
PLOTTED  ALSO  AND  ID  * 2. 

ICAL  * 1 IF  CALCOMP  PLOT  TAPE  TO  BE  CREATED.  OTHERWISE  PLOT  WILL 
BE  OUTPUT  TO  TEKTRONIX  4014  SCREEN. 

READ( 7, 1C0)NCASE,NCYCLE,ZMAX,ZMIN,XMAX,XMIN,IE,IP, IF, ID, IR, ICAL 
'.TITLE 

IF ( NCASE .EQ. 99991  GO  TO  900 

1FUFIRST  .EQ.  0)CALL  PLOTS ( I BUF  , 1000 ,22 ) 

I F ( ICAL  .EQ.  0 .OR.  IFIRST  .EQ.  1)  GO  TO  153 
CALL  PLOTIO. ,-10. ,-3) 

CALL  PLOTIO. ,1. 25,-3) 

IFIRST  • 1 
SYMBI ll-SYMI 1) 

SYMB(2)-SYMI2) 

SYMBI 3)-SYM( 3) 


PAOB  BLANC-MDt  FUND) 


■r 
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XS IZE( 1 ) *XHIN 
XSIZE(2)«XMAX 
ZSIZE(  D-ZMIN 
ZSIZEl  2)  *ZMAX 
DELTX  - ABStXMAX  - XHIN) 

DELTZ  - ABSUMAX  - ZHIN) 

I F ( I CAL  .EQ.  1 )GO  TO  154 
YRANGE  ■ DELTZ/ 11. 

XRANGE  * DEL  TX/14. 

GO  TO  157 

154  YRANGE  - DELTZ/8.5  > 

XRANGE  - DELTX/11. 

157  USE  » AMAXKYRANGE, XRANGE) 

YRANGE  » OELTZ/USE  i 

XRANGE  « OELTX/USE 

CALL  SCALE(XSIZE,  XRANGE » 2 * 1 ) 

CALL  SCALE(ZSIZE, YRANGE » 2 • 1 ) 

XP(5)«XSIZE(3> 

XP ( 6 ) *XSI ZE ( 4 ) 

ZP( 5)»ZSIZE( 3) 

ZP(6)»ZSIZE<4) 

READ  INITIAL  DATA  FROM  TAPE  3 

155  CONTINUE 

READ  (31  ICASE,  ICYCLE*  NNODE,  NELE  * TINE 
I F ( ICVCLE.GE.O)  GO  TO  200 
WRITE  (6,101) 

GO  TO  900 

200  IPLOT  « 1 

I F ( NCASE  .LT.  0)G0  TO  201 

IF< ICYCLE. £0. NCYCLE. AND. ICASE. EQ.NCASE)  GO  TO  201 
IPLOT  ■ 0 

201  WRITE  (6,102)  ICASE,  ICYCLE,  IPLOT 
READ  NODE  DATA  FRON  TAPE  3 


I F ( ID  .NE.  2 ) GO  TO  203 
DO  170  J * 1 ,NNODE 
READ( 3 ) I » X ( I ) , Z ( I ) 

X( 1+2000)  - -X(I) 

Z( 1+2000)  » Z ( I ) 

WR ITE ( 6, 103)1, X ( I ) , Z ( I I 
170  CONTINUE 
GO  TO  204 

203  READ  (3)  (NODE(I),  I-1,NN0DE) 

DO  202  J«1,NN0DE 

I « NODE ( J ) 

READ  (3)  X(I),  Y( I ) , Z(I),  IF  I X ( I ) 

I F ( IPLOT  .NE.  1 )G0  TO  202 
IF(Y(  D.NE.O.)  GO  TO  202 

WRITE  (6,103)  I,  X(I),  Y(I),  Z(I),  IFIX(I) 
202  CONTINUE 

204  IF ( IPLOT. NE.l)  GO  TO  300 
WRITE  (6,104) 


I 
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c 

C SET  UP  PLOT  AXES  IF  IPLOT  ■ l 

C ORAM  10.0  INCH  AXIS  SCALED  FOR  Z AXIS  (X  AXIS  THE  SAME  SCALE) 

C 

SUBZ3  « Z SIZE ( 3 ) 

SUBZ4  « ZS1ZE (4) 

SUBX3  * XSIZE ( 3) 

SUBX4  * XSIZE (4) 

CALL  AXIS! .5,.5,YTITLE,1 ,YRANGE ,90.0,SUBZ3,SUBZ4) 

CALL  AXIS(.5,.5,XTITLE,-1,XRANGE,0.0,SUBX3,SUBX4) 

CYCLE-NCYCLE 

CALL  TOP  (TIME, CYCLE, TITLE) 

CALL  PLOT( 5,  .5,-3) 

V C 

C READ  ELEMENT  DATA  FROM  TAPE  3 

C 

300  CONTINUE 

I F ( ID  .NE.  2)60  TO  399 
C 

C****<  TAPE  CREATED  BY  EPIC-2 
C 

DO  500  I =1 ,NELE 

READ  (3)  NEL,  Nl,  N2,  N3 , ICHECK.DI 
ISM  « 0 

IF( IPLOT. NE.l)  GO  TO  500 
IF ( ICHECK.EQ.3)  GO  TO  500 
350  CONTINUE 

C WRITE (6, 105)NEL,N1,N2,N3,ICHECK 

XP( I ) =X ( N 1 ) 

XP(2)«X(N2) 

XP( 3)«X(N3) 

XP(4)«X(N1) 

ZP( 1 )*Z(N1) 

ZP(2)-Z(N2) 

ZP ( 3 ) *Z ( N3) 

ZP ( 4 ) *Z ( N 1 ) 

C 

C CHECK  IF  TRIANGLE  VIOLATES  PLOT  BOUNDARIES  (XMIN-XMAX,  ZMIN-ZMAX) 
C 

IF ( X ( N I ) .GT. XMAX. OR. X(N2).GT. XMAX. OR. X ( N3 ) . GT.XMAX .OR . 

* Z(N1).GT.ZMAX.0R.Z(N2).GT.ZMAX.0R.Z(N3).GT.ZMAX.0R. 

* X(N1).LT.XMIN.0R.X(N2).LT.XMIN.0R.X(N3).LT.XMIN.0R. 

* Z(N1).LT.ZMIN.0R.Z(N2).LT.ZMIN.0R.Z(N3).LT.ZMIN)  GO  TO  370 
CALL  LINE(XP,ZP,4,I,0,0) 

IF ( I E .EQ.O. AND. I CHECK. EQ . 0)  GO  TO  370 
IF( IP.EQ . 0. AND. ICHECK.EQ.l)  GO  TO  370 
IF( IF. EQ.O. AND. ICHECK.EQ.2)  GO  TO  370 
C 

C WRITE  MATERIAL  SYMBOL  AT  CENTER  OF  TRIANGLE  IF  SPECIFIED 

C 

XS*( .333333*(XP(1)+XP(2)+XP(3))  - ( XS I ZE ( 3 >>) / XS I ZE < 4 ) - .02 
ZS«( .333333*(ZP(i)+ZP(2)+ZP<3))  - (ZSIZEI3) ) )/ZSIZE(4)  - .04 
ICHEC*ICHECK+ 1 
SUBSYM  > SYMB(ICHEC) 

CALL  SYMBOL! XS,ZS,. 07, SUBSYM, 0.0,1) 
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370  IF( ISM  .EQ.  1 .OR.  IR  .NE.  1>G0  TO  500 
ISW  - 1 

Nl  * Nl  ♦ 2000 
N2  * N2  ♦ 2000 
N3  * N3  ♦ 2000 
GO  TO  350 
500  CONTINUE 
GO  TO  600 
C 

C***»  TAPE  CREATED  BY  EPIC-3 

c ) 
399  DO  400  I-l.NELE 

READ  (31  NEL.  Nit  N2t  N3.  N4,  ICHECK 

IF  1 1 PLOT  .NE  • 1 ) GO  TO  400  .1 

IF( ICHECK. EO. 3)  GO  TO  400 
C 

C THREE  POINTS  ( NP1 tNP2 ,NP3)  MAKE  A PLANE 
C 

410  CONTINUE 
NP1-N1 
NP2*N2 
NP  3*N3 
IBACK-1 
GO  TO  460 
420  CONTINUE 
NP 1*N 1 
NP2*N2 
NP3-N4 
I BACK -2 
GO  TO  460 
430  CONTINUE 
NP 1 *N 1 
NP2*N3 
NP  3*N4 
IBACK-3 
GO  TO  460 
440  CONTINUE 
NP 1*N2 
NP  2*N3 
NP  3*N4 
IBACK-4 

460  IF( Y( NP1 I .EQ. 0. ANO. Y(NP2 ) .EQ.O. AND. Y(NP3) . EQ.O I GO  TO  470 
GO  TO  <420t430t440.450)tIBACt' 

450  CONTINUE 
C 

C NONE  OF  THE  PLANES  ARE  ON  THE  V-0.0  AXIS 
C 

GO  TO  400 

C 470  WRITE  (6.105)  NEL.  Nl.  N2 . N3 , N4.  ICHECK 


470  CONTINUE 
C 

C PLOT  A TRIANGULAR  TET  FACE  AT  Y « 0. 
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CHECK  IF  TRIANGLE  VIOLATES  PLOT  BOUNDARIES  IXMIN-XMAX,  ZMIN-ZMAX) 

IF(X(NPI).GT. XMA  X. OR • X ( NP2 ) • GT. XMAX. OR. X ( NP3 ) . GT .XHAX.OR. 

* Z(NP1).GT.ZMAX.0R.Z(NP21.GT.ZMAX.0R.ZINP3).GT.ZMAX.0R. 

* X(NP1I.LT.XMIN.0R.X(NP2).LT.XMIN.0R.X<NP3).LT.XMIN.0R. 

* Z(NPl).LT.ZMIN.0R.Z(K'r"2}-LT.ZNIf'-0R.Z(  NP3 ) . LT .ZMIN ) GO  TO  400 
CALL  LINE(XP,ZP«4*1,0*0) 

IF( IE.EQ.O.AND.ICHECK.EQ.O)  GO  TO  400 
I F ( IP.EQ.O.AND. ICHECK.EQ. I ) GO  TO  400 
IF( IF.EQ.O.AND. ICHECK.EQ.2)  GO  TO  400 

WRITE  MATERIAL  SYMBOL  AT  CENTER  OF  TRIANGLE  IF  SPECIFIED 


XS*< .333333*<XP(1)+XP(2)*XP(3II  - < XSI ZE (3 )!I/XS IZ E( 4 » - .02 
ZS«( .333333*(ZP(I)*ZP(2)*ZP(3))  - I LSI ZE (3 ))) /ZS IZE(4 ) - .04 
ICHEC*ICHECK+1 
SUBSYM  * SYMB ( ICHEC  > 

CALL  SYMBOL (XS,ZS,. 07, SUBSYM ,0.0,1) 

400  CONTINUE 

600  I F ( I PLOT. EQ . 0 > GO  TO  155 

GO  TO  END  OF  PLOT  AND  ESTABLISH  NEW  ORIGIN  2.0  INCHES  OVER 

I F ( ICAL  .EQ.  1 ) GO  TO  700 
CALL  PL0T(0. ,0. ,999) 

CALL  HDCOPY 

CALL  TINPUT(NHDCPY) 

CALL  PLOT!-. 5,-. 5,-3) 

GO  TO  150 

700  OVER  * OVER  ♦ 13. 

CALL  PL0T(0VER,1.25,-3) 

CALL  PLOT!-. 5, -.5,-3) 

GO  TO  150 
900  REWIND  3 

CALL  PLOT ( 0. , 0. , 999) 

STOP 
ENO 


I 


A 


c 


c 


SUBROUTINE  TOPI  TIME .CYCLE , TI TLE ) 

PLOT  TITLE  AND  HEADER  INFORMATION 

DIMENSION  TITLEI3) 

CALL  SYMBOL (2. 5*8... 14, TITLE *0.0*18) 
CALL  SYMBOL I6.0,8.*.14,6HTIME  -,0.0,6) 
CALL  NUMBER  I 999. *999. *.14, TIME *0.0 *8) 
RETURN 
END 


n r,  n c. 


r 
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APPENDIX  C 

REVISED  GPLOT  IN  EPIC-2 

SUMwUuTINt  uFLOT  (NLYCLE . I IME ) 
c: 

THIS  nOUTlNt  WHITES  DATA  UN  TmHE  3 HUH  CALCUMH  UK  It**  A01M  PLUlS 
HEVlStO  VtWSION  StPltMbtH.  1 Vif 
CUMmU.'j 


NNUDt  . 

NPNUuE . 

NELt. 

NHELt . 

ENh  W(.  y . 

o 

OT  1 . 

uTMaa , 

UTMIN, 

tmax. 

NUEH, 

ii 

PMmS. 

T MAS » 

NSLIUE. 

NC ASE  . 

SSF 

cuhmun 

a 

i)En(  10)  « 

SPn< 10) » 

E ( 10)  * 

V ( 10) » 

V I S < 10)  . 

o 

h Y ( 10)  . 

E U ( 1 0 ) * 

EU(  10)  , 

b«M ( 10)  . 

HATt ( 10)  . 

O 

CP  1 ( 1 U ) « 

CP2< 10)  ♦ 

CT1 ( 10) ♦ 

CT2< lo)  . 

tM  ( 1 0 ) * 

o 

Vb  < 10)  . 

EHS  ( 10)  . 

C ( 10)  . 

U(  10)  . 

s ( l o ) . 

o 

U<  10)  . 

tj  1 ( 1 0 ) * 

U2 1 10)  ♦ 

TEMPI ( 10)  . 

t Y Y ( 1 0 ) . 

o 

t U ( 1 0 ) » 

ELAM ( lo) , 

b< 10)  . 

A ( 10)  . 

UETVEL ( 10)  . 

o 

EusTi ( lo) « 

tuST2 ( 10)  , 

EUSTJ( 10) . 

EUST4 ( 10)  . 

E«jST5  ( 1 o ) . 

* 

EQsTto(lO) ♦ 

EUSl M 10) ♦ 

EQbTd  ( 10)  ♦ 

EusTy ( 10) * 

tuSTlo(lO) . 

« 

UAMma ( 1 0 ) . 

Et (10) 

COMMON 

« 

H ( 1000) . 

L ( 1000)  . 

HI ( 1000) * 

L l (looo)  , 

HuOT ( lDOo)  . 

o 

2UuT ( iOOO) . 

EH  ( 1000) « 

E / (1000)  . 

NOUt ( 1000  ) « 

amass ( looo) . 

o 

PMaSS ( 1000) 

IH  l*  ( 1000) 

commun 

i> 

NEL  ( (?UUU  ) » 

NOUtl (2000) 

»N0Ut2 ( 2000 ) 

.NUut  3 12000) 

. ICntCN (200u)  . 

o 

VUI.  (2000)  . 

DVUL (2000) • 

AHL A (2000) . 

MA I (2000)  . 

A I ( 2000 ) . 

i> 

t S(2000) . 

Ew (2000) . 

UArttA  (2000) 

» Du VOL (2000) 

.TSTAHT (2000) , 

HUxNt  h { ^000  ) 

CUmmQn 

IMND(S)  . 

NMNUUt (5)  . 

NSi.uUE  ( S ) » 

NStLE (b) . 

I END 1 <S)  . 

IE NU2  <S)  . 

I MN ( ^ . 50 ) . 

ISMb.SO)  ♦ 

ISt (5.30) » 

HM (50)  . 

« 

/M(bO) 

C 

lit>  FOWMAT  </////» IOX.ShCYCLE. IS. JOH  U A f A SAVtO  Ui>4  \ AHc  1 E UH  HLUl) 

C 

«HITt  (3)  NCASt.  NdrCLt.  NNODt.  NELt.HMt 
C 

DO  100  J=l.NNOUE 
I = NOUt(J) 

loO  w*I  It  (3)  1.  H(l)  » Z(i) 

C 

DU  20o  J = 1 » NELt 
I = NtL(J) 

2O0  WHITE  ( 3 > I.NOutl  < I ) ,NOUE2(I)  .NUUt3(l)  . ICHtLMi),  UOVOL  ( I ) 

C 

wHlTt  (6.1l3>  NCYCLE 
C 

Wt TUhn 
END 


4S 
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APPENDIX  D.  LISTING  OP  ISOMETRIC  PLOTTING  PROGRAM 

COMDIM*  PROC 
COMMON 

*ALPHA»NPNTS» I MAT* I START* I END 

*, SYM( 3) * SYMB ( 31 • XSI ZE (4) ,ZSIZE(5) ,XP(6) »ZP(6 ) ,NP!4) 

*, NCYCLE, NCASE 

*, TRAN SF< 4000, 4),  INTEL! 4000) 

*, NODL AY, NONOOE , I STARG , LNPROJ# I R , I TRI AN( 3200 , 2 ) 

*, ITARGT , INtLIMIT 

*«ARRAY(4,4),ARRINT(4,4) ,ARRCHGC4,4I , I ROT AT 
* *, MASK! 20000) 

END 

( 


INCLUDE  C0M0IM,LI ST 

DIMENSION  TITLE (4) ,LAYN0D(5) ,N0DE(2000) ,XYZ(4000,4) 

*t IBUF(IOOO) 

EQUIVALENCE  (MASK, NODE) » I MASK (2001) ,XYZ(1, l) ) 

DATA  LAYNOD/6, 15,24,37,50/ 

OATA  IFIRST/O/ 

C 

C ♦♦♦♦♦FORMAT  STATEMENTS***** 

C 

110  FORMAT  ( 215, 4F10. 0,11, 9X,3A6,A2) 

120  FORMAT  ( 1H1 , 5X, 7HCASE  -,I5,/,5X,7HCYCLE  «, 15,/ ,5X,7HIPL0T  »,I5) 
130  FORMAT  (1615) 

140  FORMAT  (•  INDEX  - *,110) 

180  FORMAT  (•  ITRIAN  NEEDS  TO  BE  REDI MENSIONEO  LARGER') 

C****  LIMIT  IS  DIMENSION  OF  I TRI AN*************** 

C**** 

LIMIT  « 3200 

C**** 

c 

C INITIALIZE  PLOT  AND  ESTABLISH  ORIGIN  AT  LOWER  LEFT 

C 

REWIND  3 
REWIND  7 
C 

C READ  INPUT  DATA  FOR  PLOT 
C 

C****  IOUTPT  « 1 IF  CALCOMP  PLOT  TAPE  TO  BE  GENERATED. 

C****  * 2 IF  PLOT  TO  BE  OUTPUT  TO  TEK  4014  SCREEN  IMMED. 

C****  « 3 IF  PLOT  TAPE  TO  BE  GENERATED,  LATER 

C**^*  RETRIEVED  AND  PLOTTED  ON  TEK  4014  SCREEN. 

C****  NRING  * NO.  RINGS/PROJECTILE 

C**+*  NONODE  > OR  * HIGHEST  NODE  NUMBER  FOR  TARGET 

C**»*  ISTARG  « 1ST  NODE  NUMBER  ON  TARGET 

C****  LNPROJ  « HIGHEST  PROJECTILE  NODE  NUMBER 

C**»*  IR  * 1 IF  YOU  WANT  Y REFLECTION  OF  PROJECTILE  AND  TARGET 
C**»*  CREATED.  OTHERWISE  IR  ■ 0 

C^***  ITARGT  * 1 IF  THIS  PLOT  INCLUDES  A TARGETdE.  HIDOEN 
C****  LINE  ALGORITHM  MUST  BE  USED). 

(;♦♦♦♦  NX  = NUMBER  OF  TARGET  ELEMENTS  ALONG  X AXIS 

C****  NY  * NUMBER  OF  TARGET  ELEMENTS  ALONG  Y AXIS 

C****  NZ  * NUMBER  OF  TARGET  ELEMENTS  ALONG  Z AXIS 

183  READ  (7,110)  NCASE , NCYCLE ,ZMAX, ZMI N,XMAX,XMI N, IOUTPT , T ITLE 
I F INCASE  .EQ.  9999)G0  TO  610 

READ  (7,130)  NRING, NONODE, LNPROJ, ISTARG, IR, ITARGT, NX, NY 
*»NZ 

C***» INITIALIZE  IN  FOR  CKTRl  SUBROUTINE.  IN  * NO.  TRIANGLES 
£♦♦♦♦  ALREADY  CONSIDERED  ♦ 1 
IN  * 1 

C****  IMAT  ■ 1 FOR  PROJECTILE 
C****  » 2 FOR  TARGET 


C****  IMAT  ■ 
O***  * 

IMAT  * 
C****  NODLAY 
NODLAY 


NUMBER  NODES/LAYER  OF  PROJECTILE 
LAYNOD(NRING) 


DO  187  I > 1,4000 
INTEL! I ) ■ 0 


1 


46 


o o o 


? 


00  186  L > 1*4 

186  TRANSFC I . L ) « 0. 

187  CONTINUE 

1M0D *LAYNOD (NR1NG  - 1) 

JMOO-IMOD  ♦NONODE 
KMOD«NOOLAY*NONODE 

C****  THE  FOLLOWING  ASSUMES  2 LE  NRI NG  LE  5 
ICKl  « ( NRING  - 2)  * 2 ♦ 6 
ICK2  * ICKl  ♦ 2 

C****  SEi  IN01  AND  IN02  FOR  CKPTS  SO  THAT  NODES  1 THRU  LNPROJ  ARE 
C**»*  CHECKED  AGAINST  EXTERIOR  TRIANGLES  OF  PROJECTILE 
IND1  - 1 
IND2  * LNPROJ 
IRELEM«2*N0N0DE 

IREF  « NONODE  * 10000  ♦ NONOOE 
IF  (IR.EQ.O)  IRELEM-NONODE 
185  DO  190  I »1 * IRELEM 

XYZ( I,4)»l 
190  CONTINUE 

C 

SYMB ( 1 1 *SYM(  1 ) 

SYMB ( 2 ) *SYM ( 2 ) 

SYMBI 3 1 *SYM ( 3 I 

C****  BASED  ON  11  X 14  TEKTRONIX  4014  SCREEN  DETERMINE  MAXIMUM 
C«***  SIZE  PLOT  CAN  BE  TO  FIT  ON  SCREEN 
XS IZE( 1 l =XMIN 
XSIZEI 21 *XMAX 
ZSIZEI 1)«ZMIN 
ZSIZEC21 *ZMAX 
DELTX«ABS( XMAX-XMIN) 

DELTZ*ABS(ZMAX-ZMIN» 

IFCIOUTPT  .EQ.  1 ) GO  TO  195 
YRANGE*DELTZ/11. 

XRANGE*DELTX/14. 

GO  TO  197 

195  YRANGE  * OELTZ/8.5 
XRANGE  « DELTX/11. 

197  USE«AMAX1(YRANGE, XRANGE) 

YRANGE*DELTZ/USE 
XRANGE *DELTX/USE 
CALL  SCALE  ( XSI ZE * XRANGE ,2 *1 ) 

CALL  SCALE  ( ZSI ZE .YRANGE *2 . 1 t 
XP(5)«XSIZE( 3) 

XP ( 6 ) *XS I ZE ( 4) 

ZP( 5)*ZSIZE( 3) 

ZP(6)-ZSIZE(4) 

READ  INITIAL  DATA  FROM  TAPE  3 

READ  (3)  ICASE * IC YCLE .NNODE  *NELE  * TI ME 
IF  ( ICYCLE.GE.OI  GO  TO  210 
WRITE  (6*101) 

STOP 
IPLOT-1 

IF  ( ICYCLE.EQ.NCYCLE. AND. ICASE. EQ.NCASE)  GO  TO  220 


200 


210 
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IPLOT-O 

WRITE  (6.1201  ICASE « (CYCLE* I PLOT 
READ  NODE  DATA  FROM  TAPE  3 


READ  (31  (NODE(I»,I«l,NNOOE) 

DO  230  J«1 .NNOOE 
I -NODE! J I 

READ  (3)  XYZ ( I * 1 1 , XYZ (I«2I,XYZ(I,3),IFIX 
230  CONTINUE 

IF  (IPLOT.EQ.OI  GO  TO  300 

IF  ( IR.EQ.OI  GO  TO  250 

C****  CREATE  MIRROR  IMAGE  OF  PENETRATOR 
DO  240  J*l. NNOOE 
I-NODE(J) 

XYZ(UNONOOE»U«XYZ(I  *1) 

XYZ ( I+NONOOE, 2 1 »-XYZ(  I .2) 

XYZ( I+N0N00E»3)»XYZ(I  .3) 

240  CONTINUE 

C****  READ  IN  SCALING . TRAN SLAT  I ON .ROT  ATI ONAL  PARAMETERS 
250  CALL  TRANS 
C****  TRANSFORM  COORDINATES 
C 

DO  280  I «1 » IRELEM 
DO  270  L*1 ,4 
DO  260  K*1 ,4 

TRANSF(I ,L)»XYZ(I ,K) 4ARRI NT (K, L) +TRANSF ( I , L ) 

260  CONTINUE 

270  CONTINUE 

280  CONTINUE 

283  ALPHA  * 3.14159/2. 

C**«*  INITIALIZE  MASK  SO  THAT  TARGET  WILL  NOT  BE  DRAWN  WHERE  IT 
C**«*  COINCIDES  WITH  PROJECTILE. 

C 

285  DO  290  K«l, 20000 

MASK (K)*- 10000 
290  CONTINUE 

CYCLE-NCYCLE 

I F ( IOUTPT  .NE.  1 .OR.  IFIRST  .EQ.  01 
♦CALL  PLOTS( IBUF, 1000,22) 

IF( IOUTPT  .NE.  1 .OR.  IFIRST  .EQ.  1 )G0  TO  295 
IFIRST  = 1 

CALL  PL0T(0.,-10. ,-3) 

CALL  PL0T(0.»1.25,-3) 

295  IF( IOUTPT  .EQ.  3ICALL  PLOT ( 0. ,0. ,997 ) 

CALL  TOP  (TIME, CYCLE, TITLE) 

C 

C READ  ELEMENT  DATA  FROM  TAPE  3 
C 

300  DO  600  IA*1,NEL£ 

READ  (3)  NEL ,N1 ,N2 ,N3 ,N4 , ICHECK 
IF  (1PL0T.EQ.0)  GO  TO  600 

C***»  CHECK  TO  SEE  IF  FIRST  TARGET  NOOES  READ  YET,  IF  SO 
£*♦♦♦  PROCESS  PROJECTILE  NODES 

IF(N1  .GT.  LNPROJ  .AND.  IMAT  .EQ.  1 )G0  TO  315 
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C****CHECK  4 TRIANGLES  OEFINED  BY  4 POINTS  READ 
C***« IGNORE  INTERIOR  ONES. 

305  CALL  CKTRI  (N1 ,N2 ,N3,N4) 

IF( I A .NE.  NELE )G0  TO  600 
C****WHEN  LAST  SET  OF  POINTS  READ  OR.  IF  TARGET, 

C****F IRST  SET  OF  TARGET  POINTS  READ,  CONE  HERE. 

315  CONTINUE 

IF  ( IR.NE.l)  GO  TO  330 

C****PACK  ITR IAN  SO  THERE'LL  BE  ROOM  FOR  REFLECTIVE  TRIANGLES. 

CALL  PACK  (INDEX) 

WRITE  (6,140)  INDEX 

IF  ( ( INDEX* INDEX) .GT. LIMIT)  GO  TO  520 
C****IF  Y REFLECTION  OF  PROJECTILE  AND  TARGET  DESIRED  CREATE 
C****RE ELECTIVE  TR I ANGLES( EXTERIOR  ONLY). 

DO  320  I«l, INDEX 

ITR IAN ( I* INDEX, 1 )*ITRIAN(I,1)*IREF 
ITRIAN(I  ♦ INDEX, 2)  * ITRIAN(I,2)  ♦ NONODE 
320  CONTINUE 

C****CHECK  ALL  POINTS  AGAINST  EXTERIOR  TRIANGLES  TO  SEE  WHICH 
C****POINTS  ARE  HIDOEN  AND  SHOULD  NOT  BE  CONNECTED. 

330  CALL  CKPTS(INDl,I ND2 ) 

NPNTS  * 4 

DO  490  JJ*1. LIMIT 
IF  ( ITRIANI JJ,l).EQ.O)  GO  TO  490 
C****RETR I EVE  NODAL  POINTS  FROM  ITRIAN 
Ml  * I TR I AN ( J J , 1 ) / 10000 
M2  * ITR I AN( J J , 1 ) - (Ml  * 100001 
M3  > ITRI AN( JJ,2) 

C+**+ I F ONE  OF  THE  THREE  POINTS  DEFINING  A TRIANGLE  IS  HIDDEN, 

C****NO  PART  OF  TRIANGLE  IS  PLOTTED. 

IF( INTEL (Ml)  .EQ.  1 .OR.  INTEL (M2)  .EQ.  1 .OR. 

* INTEL ( M3 ) .EQ.  1)G0  TO  490 

C****SINCE  A TRIANGLE  IS  PLOTTED  WHEN  ALL  3 POINTS  DEFINING  IT  ARE  SEEN, 
C****THERE  ARE  SOME  ANGLES  AT  WHICH  AN  ENTIRE  PROJECTILE 
C****!S  PLOTTED  WHEREIN  SOME  LINES  AT  THE  FAR  END  OF  THE 
C****PROJECTILE  ARE  UNDESIRABLE. 

C**** 

IF( IMAT  .EQ.  2 )G0  TO  340 
I F < IR  .EQ.  1)G0  TO  407 
C 

C***»  PROJECTILE  ONLY  - NO  REFLECTION 
C 

GO  TO  410 
C 

€***•  HALF  AND  FULL  TARGET 

C 

C****  THERE  ARE  ALSO  SOME  ANGLES  AT  WHICH  THE  TARGET  IS 
C****  PLOTTED  WHEREIN  SOME  LINES  NEAR  CORNERS  ARE  UNDESIRABLE 
340  KT1  - Ml  - ISTARG  ♦ 1 
KT2  - M2  - ISTARG  ♦ 1 
KT 3 ■ M3  - ISTARG  ♦ 1 
NXNY  ■ NX  * NY 
NXNYM  ■ NXNY  ♦ ISTARG  - NX 
NXNYMR  > NXNYM  ♦ NONOOE 

C**«*  00  NOT  PLOT  TRIANGLE  JOINING  UPPER  LEFT  BACK  CORNER 


SI 


C****  OF  TARGET  TO  POINT  BELOW  IT  ANO  TO  ITS  RIGHT 

IF< IABSIM1-M2)  .EQ.  1 .ANO.  IABS(M1-M3)  .EQ.  NX NY 

♦ .AND.  MOD(KTl.NX)  .EQ.  1 .ANO.  MOD(KT3.NX>  .EQ.  1 

♦ .ANO.  HI  .EQ.  NXNYHI 
♦GO  TO  490 

IF( IR  .EQ.  0)G0  TO  405 
C 

FULL  TARGET 
C 

C***»  REFLECTIVE  TRIANGLES  CREATED  FOR  1ST  NX  NODES  OF  TARGET 
C****D0  NOT  SHOW  UP  AS  DUPLICATE  INTERIOR  TRIANGLES  BECAUSE 
C****THE  NODE  NUMBERS  ARE  DIFFERENT  SO  CONNECTIONS  JOINING  1ST 
C****NODE  OF  TARGET  (OR  ITS  REFLECTIONITO  NODE  BELOW  AND  TO  THE 
C****RIGHT  MUST  BE  ELIMINATED. 

I F C IABSI M1-M2 ) .EQ.  1 .AND.  IABSIM1-M3)  .EQ.  NX NY 

♦ .AND.  (Ml  .EQ.  ISTARG  .OR.  Ml  .EQ.  (ISTARG  ♦ NONODE ) ) ) GO  TO  490 
I F ( M 1 .IE.  NONODE (GO  TO  410 
I T 1 * KT1  - NONODE 
IT2  * KT2  - NONODE 
IT3  ■ KT3  - NONODE 

C ♦♦♦♦  DO  NOT  PLOT  TRIANGLE  JOINING  UPPER  RIGHT  FRONT  CORNER  OF 
C**«*  FULL  TARGET  TO  POINT  BEHIND  IT  AND  BELOW  IT(NOTABLY  FOR 
C****  POSITIVE  ROTATION  FROM  X TO  Y AXIS 

IF( IABSI Ml  - M2)  .EQ.  NX  .AND.  IABS(H2  - M3) 

♦ .EQ.  NXNY 

♦ .AND.  IT2  .EQ.  NXNY )G0  TO  490 
C****  DO  NOT  PLOT  TRIANGLE  JOINING  LOWER  RIGHT  FRONT  CORNER 
C *♦♦♦  OF  TARGET  TO  POINT  BEHIND  IT  AND  POINT  ON  ITS  LEFT 

IF( IABSIM2-M3)  .EQ.  1 .AND.  IABSIM1-M3)  .EQ.  NX  .AND. 

♦ ITl  .GT.  NXNY)GO  TO  490 
C***«  UPPER  LEFT  FRONT  CORNER 

I F ( IABSI M 1-M2 ) .EQ.  NX  .AND.  IABSIM2-M3)  .EQ.  NXNY 

♦ .AND.  M2  .EQ.  NXNYMR ) GO  TO  490 
GO  TO  410 

C 

(;♦♦♦♦  HALF  TARGET 

C 

C****  UPPER  LEFT  CORNER 

405  I F ( IABS( M1-M2)  .EQ.  NX  .AND.  IABS(M1-M3)  .EQ.  NXNY 

♦ .AND.  MOD ( KT1 » NX ) .EQ.  1 )G0  TO  490 
C+*+*  LOWER  RIGHT  CORNER 

IF ( I ABS ( M 1-M2 ) .EQ.  1 .AND.  IABSIM2-M3)  .EQ.  NX 

♦ .AND.  XT  1 .GT.  NXNY)GO  TO  490 
C ♦♦♦♦  UPPER  RIGHT  CORNER 

IF ( IABS(M1-M2)  .EQ.  NX  .AND.  IABSIM1-M3)  .EQ.  NXNY 

♦ .AND.  MOD(KTl.NX)  .EQ.  0 .AND.  M0D(KT2.NX)  .EQ.  C 

♦ .AND.  MOD( KT3» NX ) .EQ.  0) 

♦ GO  TO  490 

C****  LOWER  LEFT  CORNER 

IF ( IABS(M1-M2)  .EQ.  1 .AND.  IABS(M1-M3)  .EQ.  NX 

♦ .AND.  MOD(KTl.NX)  .EQ.  1 

♦ .AND.  KTl  .GT.  NXNY)GO  TO  490 
GO  TO  410 

C 

C**«*  PROJECTILE  ONLY  - WITH  REFLECTION 


' 


52 


I 
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407  I F ( M 1 .LE.  IM00  .AND.  (M2  .GT.  NODLAV  .OR.  M3  .GT.  NO OLAY  I ) 
♦GO  TO  490 

IFIMl  .GT.  NONODE  .AND.  Ml  .LE.  JMOO  .ANO. 

♦ (M2  .GT.  KMOD  .OR.  M3  .GT.  KMOD) I GO  TO  490 
410  NP(1)«M1 

NP(2)«M2 
NP ( 3) -M3 
NP(4)-M1 

420  DO  430  I -1,NPNTS 

J-NP(I) 

XP(I)-TRANSF( J,l) 

ZP( I)-TRANSF( J,3> 

430  CONTINUE 

C WRITE(6,4)(XP(I),ZP(I) , 1-1,6) 

C 4 FORMAT( 12F9.3) 

431  I F ( ITARGT  .EQ.  0)G0  TO  435 
CALL  PL0T30 
GO  TO  490 

C****IF  NO  TARGET,  PLOT  PROJECTILE  WITHOUT  SETTING  UP 
^♦♦♦HIDDEN  LINE  MASK  FOR  TARGET. 

435  KK  - 3 

DO  440  I - 1 , NPNTS 

AX  = (XP( I)-XP(5) )/XP(6) 

AY  = (ZP(I)-ZP(5) )/ZP(6) 

CALL  PLOT( AX, AY,KK) 

KK  3 2 
440  CONTINUE 
490  CONTINUE 

00  495  I « l, LIMIT 
495  ITR I AN( 1,1)  * 0 

I F ( ITARGT  .NE.  1 .OR.  IMAT  .EQ.  2)G0  TO  605 
C****  INITIALIZE  IN  FOR  CKTRI 
IN  = 1 
IMAT  * 2 
CALL  CKMASK 

C**«*  INITIALIZE  IND1  AND  IN02  FOR  CKPTS 
IN01  - ISTARG 
IN02  * NONOOE 
GO  TO  305 
600  CONTINUE 

IFIIPLOT  .EQ.  0)G0  TO  200 
605  I F ( IOUTPT  .EQ.  1)G0  TO  608 
CALL  PL0T(0. 0,0. 0,999) 

I F ( IOUTPT  .EQ.  3 ) GO  TO  183 

CALL  HOC  OP Y 

CALL  TINPUT(IOUM) 

GO  TO  183 

608  OVER  - OVER  ♦ 13. 

CALL  PLOT) OVER, 1. 25,-3) 

GO  TO  183 
610  REWIND  3 

CALL  PL0T(0.,0.,999) 

STOP 

520  WRITE  (6,180) 

STOP 

END 
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SUBROUTINE  TRANS 
INCLUDE  COHDIM 
30  CALL  IDENT 
200  READ! 7*1) INORD 
l FORMAT! Ill 
C****  SCALING 

I F ( I WORD  .EQ.  1)  GO  TO  300 
C**»*  TRANSLATION 

I F C IWORD  .EQ.  2)  GO  TO  400 
C**»*  ROTATION 

IF ( IWORD  .EQ.  3)  GO  TO  500 
C****  PLOT 


I F ( IWORD  .EQ.  4)  RETURN 
C«***  START  OVER 

IF( INORD  .EQ.  5)  GO  TO  30 
GO  TO  200 

300  READ(7,2)IS1,IS2.IS3 
2 FORMAT! 6141 

51  * IS1/100. 

52  * IS2/100. 

53  * IS3/100. 

CALL  SCALIN! S1,S2,S3) 

GO  TO  200 

400  READ!7,2)IT1,IT2.IT3 
T 1 = IT1/100. 

T2  * IT2/100. 

T3  * IT3/100. 

CALL  TRANSHT1»T2»T3) 

GO  TO  200 

500  READ!7,2>IRlt IR2.1THETA 
550  CALL  ROTATE! IR 1 , IR2 , 1 THE TA I 
GO  TO  200 
END 


c **** 
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SUBROUTINE  IDENT 

INITIALIZE  TRANSFORMATION  MATRICES 

INCLUDE  COMDIM 

00  10  I > 1,4 

DO  5 J > 1,4 

ARRINTI I , J)  * 0 

ARRAY ( I , J I > 0 

ARRCHGI I , J ) = 0 

CONTINUE 

DO  20  I * 1,4 

ARRAY! 1,1)  * 1 

ARRINT! I , I I * 1 

RETURN 

END 


SUBROUTINE  SCALIN! SI tS2,S3) 

€****  PERFORM  SCALING 
INCLUDE  COMOIM 
ARRAY! 1« I)  > SI 
ARRAY! 2, 21  » S2 
ARRAY ( 3 * 3 > - S3 
CALL  HULMAT 
RETURN 
ENO 

\ 


SUBROUTINE  ROTATE ( N1 ,N2 . I THETAJ 
C****  PERFORM  ROTATION  OF  MATRICES 
INCLUOE  COMDIM 
THETA  = ITHETA 

STHETA  = SIN ( THETA  * .017*532931 
CTHETA  » COS( THETA  * .017*532931 
ARRAY ( N 1 1 N1 ) ■ CTHETA 
ARRAY(N2«N2 ) « CTHETA 
ARRAY(NlfN2)  ■ STHETA 
ARRAY(N2,Nl)  - - STHETA 

CALL  MULMAT  ) 

RETURN 

END 

) 
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SUBROUTINE  MULMAT 
C+***  MULTIPLY  MATRICES 
INCLUDE  COMDIM 
DO  10  I * 1,4 
DO  5 J * 1,4 
DO  3 K > 1,4 
ARRCHG! I , J ) * 
CONTINUE 
CONTINUE 
DO  30  I = 1,4 
DO  20  J * 1,4 
ARR INTI  I , J ) ■ 
ARRAY ( I , J ) « 0 
ARRCHGI I , J ) * 0 
CONTINUE 
DO  40  I > 1,4 
ARRAY! I, I)  » 1 
RETURN 
END 


3 

5 

10 


20 

30 

40 


ARRAY! I , K)  * ARRINTIK, J)  ♦ ARRCHGII.JI 


ARRCHG ( I , J) 


SUBROUTINE  TOP ( TI HE *C YCLE .TITLE I 

PLOT  TITLE  AND  HEADER  INFORMATION 

DIMENSION  TITLEOI 

CALL  SYMBOL (2.5,8.,. 1A,TITLE, O.O, 18) 
CALL  SYMBOL (6.0,8. ,.14,6HT1ME  >,0.0,6) 
CALL  NUMBER ( 999., 999. , . 1 4 , TI ME , 0. 0 ,8 ) 
RETURN 
END 


SUBROUTINE  CKTRI  INI ,N2,N3,NA) 

C****  THIS  ROUTINE  CHECKS  FOR  INTERIOR  TRI ANGLES  I SHOULD  BE 
C****  DUPLICATED)  AND  DROPS  THEM. 

INCLUOE  COMDIM.LIST 
DIMENSION  NUM!A,2),NN!A) 

DATA  MULTI/IOOOO/ 

C 

C *****FORMAT  STATEMENTS***** 

C 

110  FORMAT  !•  ITRIAN  NEEDS  TO  BE  REDI MENSI ONEO  LARGER1) 

C 

IN  = IN 

NN! 1)«N1 

NNI2J-N2 

NN  ( 3 ) *N3 

NN ( A ) *NA 

IC0UNT*IC0UNT*1 

C****  SORT  NODAL  POINTS  SO  THAT  N1<N2<N3<NA 
M«2 

DO  130  K * 1 » 3 
DO  120  L*M,  A 

IF  (NN(L).GT.NN(K) ) GO  TO  120 
I SAVE  *NN ( K ) 

NN!K)»NN!L) 

NN ( L ) * I SAVE 
120  CONTINUE 

M*M+1 

130  CONTINUE 

C****  PACK  THE  NODAL  NUMBERS  OF  THE  THREE  POINTS  MAKING  UP  TRIANGLE 
C**»* INTO  1 WORD  IN  ITRIAN.  IF  POINTS  1000,2000,3000  MAKE  UP  TRIANGLE, 
C****ITRIAN(  1,1)  =>  10002000  AND  I TR  I AN(  1 ,2 ) « 3000 
NUM(l.l)  * NN ( 1 ) * MULTI  ♦ NN<2) 

NUM( 1,2)  * NN ( 3 ) 

NUM( 2,1)  » NUM (1,1) 

NUM( 2,2)  » NN ( A ) 

NUM( 3,1)  * NN ( 1 ) * MULTI  ♦ NN(3) 

NUM ( 3,2)  * NN ( A ) 

NUM( A, 1 1 = NN ( 2 ) * MULTI  ♦ NN(3) 

NUM( A, 2 ) » NN( A) 

C****  CHECK  THIS  TRIANGLE  AGAINST  OTHERS  ALREADY  PROCESSED.  IF  A 
C****  TRIANGLE  IS  PROCESSED  TWICE,  IT  MUST  BE  AN  INTERIOR  TRIANGLE, 
C****  SO  IT  IS  REMOVED  FROM  ITRIAN  ( I TRI AN( I , 1 ) *0 ) AND  NOT  CHECKED 
C****  AGAINST  POINTS  TO  SEE  IF  IT  HIDES  THEM. 

DO  170  K*l,A 

DO  1A0  L«1,IN 

IF( I TR I AN ( L , 1 ) .NE.  NUM(K,1))G0  TO  1A0 
I F ( I TR I AN( L , 2 ) .EQ.  NUM(K,2))G0  TO  150 
1A0  CONTINUE 

ITRIAN! IN, 1)>NUM(K,1) 

ITRIAN! IN, 2)  « NUM!K,2> 

IN-IN+1 

IF  ! IN. GT. LIMIT)  GO  TO  160 
GO  TO  170 

150  ITRIANtL, 1 (*0 

GO  TO  170 


61 


c**»* 

160 


1 70 


ITRIAN  IS  DIMENSIONED  TO  LIMIT.  IF  AN  OVERFLOW  Ofr.URS,  GO 
BACK  AND  FILL  IN  SPACES  FROM  WHICH  TRIANGLES  WERr  REMOVED. 
CALL  PACK  (INDEX) 

IN=INOEX 

IF  ( IN. GT. LIMIT)  WRITE  (6,1101 
CONTINUE 
RETURN 
END 


SUBROUTINE  PACK! INDEX) 

C****  THIS  ROUTINE  PACKS  I TR I AN.  BLANK  SPOTS  WERE  CREATED  WHEN  DUPLICATE 
C****  TRIANGLES  WERE  REMOVED. 

INCLUDE  COMDIM.LIST 
M « LIMIT 
DO  10  J * 1, LIMIT 
!F< ITRIANI J,l)  .NE.  OJGO  TO  10 
00  5 I * M, 1,-1 
IF( I TR I AN (1,1)  .NE.  0)G0  TO  8 
5 CONTINUE 
GO  TO  20 

8 I F ( J .GT.  I )G0  TO  20 

ITRIANIJ.l)  = I TR I AN ( I ,1 ) 

ITR I AN ( J , 2 ) = I TR I AN  C I ,21 
I TR I AN (1,1)  « 0 
M « I - 1 

IFU  ,GE.  M )G0  TO  15 
10  CONTINUE 
15  INDEX  - J 
GO  TO  25 
20  INDEX  * J - 1 
25  WR ITE< 6, 1>( ITRIANI 1,1), I *1, INDEX) 

1 FORMAT! 2120) 

RETURN 

END 
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SUBROUTINE  CKPTS ( IND1 , I N02 ) 

C**»*  THIS  ROUTINE  CHECKS  ALL  POINTS  DEFINING  EXTERIOR  TRIANGLES 
C****  AGAINST  ALL  EXTERIOR  TRIANGLES  TO  SEE  WHICH  POINTS  ARE 
C****  HIDDEN  and  should  not  be  connected. 

INCLUDE  COMDIM,LI ST 
C 

C *****FORMAT  statements***** 

c 

120  FORMAT  (6F10.3,I2) 

C 

DO  480  JJ*1, LIMIT 
IF  ( ITRIANI JJtlI.EQ.O)  GO  TO  480 
MI  = ITRIANI  JJ, 11/10000 
M2  = I TR I AN ( J J « 1 ) - (Ml  * 10000) 

M3  = ITR I AN ( J J • 2 ) 

C****IF  ALL  3 POINTS  OF  TRIANGLE  HAVE  SAME  X OR  Z COOROI NAT ES , DO 
C****NOT  CONSIDER  THIS  TRIANGLE. 

I F ( TRANSF ( M 1 t 1 ) .EQ.  TRANSF ( M2 • 1 ) .AND. 

* TRANSF (Ml* 1 ) .EQ.  TRANSF ( M3 , 1 )) GO  TO  480 
I F ( TRANSF ( Ml » 3 ) .EQ.  TRANSF ( M2  1 3)  .AND. 

* TRANSF ( Ml « 3 ) .EQ.  TRANSF ( M3 , 3 ) ) GO  TO  480 
C****!F  2 POINTS  OF  TRIANGLE  HAVE  SAME  Z AND  X COORDINATES, 

C****DO  NOT  CONSIDER  THIS  TRIANGLE. 

IF  ( < TRANSF ( Ml , 3). EQ. TRANSF (M2, 3). AND. TRANSF (Ml, II 

1 .EQ. TRANSF (M2, 1) I.OR. ( TRANSF ( Ml , 3 ) . EQ.TR ANSF ( M3, 3 1 

2 . AND. TRANSF (Ml, ll.EQ. TRANSF! M3, 1)). OR. ( TRANSF ( M2, 3 ) 

3 .EQ.TRANSF(M3,3).AND.TRANSF(M2,l).EQ.TRANSF(M3,im  GO  TO  480 
C****DETERMINE  1ST  POINT  OF  TRIANGLE 

C****  THIS  IS  ARBITRARILY  ONE  WITH  LOWEST  Z VALUE. 

IF  (TRANSF(M1,3I-TRANSF(M2,3I)  150,140,130 
130  IF  (TRANSF(M2,3)-TRANSF(M3,3>)  190,200,160 

140  IF  (TRANSF(M1,3)-TRANSF(M3,3)>  170,480,160 

150  IF  (TR AN SF (Ml, 3)-TRANSF( M3, 31  I 180,210,160 

160  IPT1*M3 

IPT2*M2 
IPT3*M1 
GO  TO  220 

170  IF  (TRANSF(M1,1)-TRANSF(M2,1>)  180,480,190 

180  IPT1*M1 

IPT2*M2 
IPT3*M3 
GO  TO  220 
190  IPT1*M2 

I PT2*M3 
IPT3-MI 
GO  TO  220 

200  IF  ( TRANSF ( M2 » l l-TRANSF ( M3  » l ) ) 190,480,160 

210  IF  ( TRANSF(M1, l l-TRANSF (M3 ,1 ) I 180,480,160 

C****  DETERMINE  SLOPES  BETWEEN  POINT  1 AND  OTHER  2 POINTS 
C****IN  ORDER  TO  DETERMINE  POINTS  2 AND  3 SO  THAT  1-3 
C****DEFINES  TRIANGLE  IN  COUNTER  CLOCKWISE  ORDER. 

C****DETERM INE  SLOPE  BETWEEN  POINTS  1 AND  2. 

220  IF  (TRANSF(IPT1,1)-TRANSF(IPT2,1)I  240,230,240 

230  AM  1*99999. 

GO  TO  250 


I 
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240  AM1«( TRANSF ( IPT1,3)-TRANSF ( I PT2  ,3))/(TRANSF( IPTltll- 

1 TRANSF! IPT2, 111 

IF  ( ABSIAMD.LT. .0001 ) AMI-0.0 
C**»* DETER MINE  SLOPE  BETWEEN  POINTS  1 ANO  3. 

250  IF  ( TRANSF ( I PT1 1 1 1 -TRANSF ( IPT3 « D ) 270*260*270 
260  AM3-99999. 

GO  TO  280 

270  AM3«( TRANSF I IPT1 ,31-TRANSF ( I PT3,3> )/ (TRANSF!  IPT1.1I- 

1 TRANSF!  IPT3*  ID 

IF  (ABS( AM3I.LT.. 0001)  AM3-0.0 

C****IF  BOTH  SLOPES  ARE  POSITIVE  OR  BOTH  ARE  NEGAT IVE*T AKE  POINT 
C****WHOSE  SLOPE  WITH  POINT  1 IS  SMALLER. 

280  IF  ( (AMI. GE.O.. ANO. AM3.GE.0.). OR. (AMI. LT. 0.0. AND. AM3.LT. 0.) ) 

1 GO  TO  290 

C**** I F ONE  SLOPE  IS  POSITIVE*THE  OTHER  NEGATIVE  * POI NT  2 IS  ONE 
C****WHOSE  SLOPE  WITH  P0INT1  IS  POSITIVE. 

IF  (AMI. GE.O.)  GO  TO  310 
GO  TO  300 

290  IF  (AM1.LT.AM3)  GO  TO  310 

300  $AVE*AM3 

I SAVE  3 IP  T3 
AM3-AM1 
AMI-SAVE 
I PT3* I PT2 
IPT2* I SAVE 
310  CONTINUE 

C****DETERMINE  SLOPE  BETWEEN  POINTS  2 AND  3. 

IF  ( TR  AN  SF(IPT2,1)- TRANSF  (IP  T3,D)  330,320,330 
320  AM2-99999. 

GO  TO  340 

330  AM2-( TRANSF (IPT2,3)-TRANSF(I PT3 »3))/(TRANSF(IPT2»l)“ 

1 TRANSF! IPT3, l) ) 

IF  ( ABS( AM2) .LT..0001 ) AM2-0.0 
C**»*  DETERMINE  WHICH  POINTS  ARE  HIDDEN  BY  THIS  TRIANGLE. 

C**»*  IF  HIDDEN,  INTEL ( I ) - 1 

340  LL-IND1 

MM-IND2 


DO  470  KK *1 *2 

DO  460  DLL, MM 

C****IF  POINT  ALREADY  HIDDEN  BY  A TRIANGLE  DONT  TEST  IT  AGAIN. 

IF  (INTEL(I).EQ.l)  GO  TO  460 

C****IF  Y( POINT ) <0R*  Y1,Y2,Y3,  IT  CANT  BE  HIDDEN  BY  TRIANGLE 

IF  ( TRANSF ( I ,2 ) . LE. TRANSF ( I PT1 ,2  I . AND.TRANSF ( 1,2) 

1 .LE. TRANSF (I PT2 ,2 ). AND. TRANSF (1,2). LE. TRANSF! IPT3.2)  ) 

2 GO  TO  460 

C****IF  X( POINT )> OR*  X1,X2,X3,  IT  CANT  BE  HIDOEN  BY  TRIANGLE 

IF  ( TRANSF ( I , 1 ).GE. TRANSF ( I PT1,1). AND. TRANSF! 1,1) 

1 . GE. TRANSF ( I PT2 , 1). AND. TRANSF ( 1, 1).GE. TRANS F< IPT3, 1) ) 

2 GO  TO  460 

C****IF  X( POINT<OR  = X 1 , X2 , X3 , IT  CANT  BE  HIDDEN  BY  TRIANGLE 

IF  (TRANSF (I ,l).LE. TRANSF (I PT1,1). AND. TRANSF( 1,1) 

1 .LE .TRANSF ( I PT2 , 1 ) .AND.TRANSF (1,1). LE. TRANSF ( I PT3, 1 > ) 

2 GO  TO  460 

C****IF  Z ( POINT )>0R*  Z1,Z2,Z3,  IT  CANT  BE  HIDDEN  BY  TRIANGLE 

IF  ( TRANSF ( I, 3). GE. TRANSF! I PT1, 3 ). AND. TRANSF (I, 3) 
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1 

2 


C****IF  Z I POINT 


C****IF  POINT  = 


1 

2 


C****IF  POINT 


C****IF  POINT  * 


1 

2 


C****  DETERMINE 

350 


360 


1 


370 

380 


390 


400 


IF 

IF 

IF 

IF 


410 

420 


430 

440 


C***«  DETERMINE 


1 

2 
3 


1 


• GE  . TRANSF ( I PT2  *3) .AND. TRANSF I I ,3 ) . GE.T RANSF ( IPT3,  3 ) » 

GO  TO  460 

KOR ■ Z1,Z2,Z3,  IT  CANT  BE  HIDDEN  BY  TRIANGLE 
IF  (TRANSFII ,3). LE. TRANSF II PT1,3). AND. TRANSF! 1,3) 

•LE .TRANSF ( I PT2 , 3) .AND. TRANSF (1,3). LE • TRANS F ( IPT3, 3 ) I 
GO  TO  460 

P0INT1  OF  TRIANGLE.IT  CANT  BE  HIDDEN  BY  TRIANGLE 
IF  ( TRANSF ( I, 1 I. EQ. TRANSF! I PT1,1). AND. TRANSF! 1,2) 

. EQ. TRANSF ( I PT1, 2). AND. TRANSF (I, 3 ).EQ. TRANSF! IPT1.3) ) 

GO  TO  460 

P0INT2  OF  TRIANGLE,  IT  CANT  BE  HIDDEN  BY  TRIANGLE 
IF  ( TRANSF ( I , I ) .EQ. TRANSF ! I PT2 , 1 ) . AND. T RANSF! 1,2) 

. EQ. TRANSF ! I PT2, 2 ). AND. TRANSF! I, 3). EQ. TRANSF! IPT2.3)  ) 

GO  TO  460 

P0INT3  OF  TRIANGLE,  IT  CANT  BE  HIDDEN  BY  TRIANGLE 
IF  ( TR AN SF  ( I , I ) . EQ.  TRANSF  < I PT3  , l ) . AND . T R ANS F 1 1 , 2 ) 

.EQ. TRANSF! I PT3,2) .AND. TRANSF! 1,31 .EQ. TRANSF ! IPT3,  3) ) 

GO  TO  460 
I SL0SW=0 

SLOPES  BETWEEN  POINT  AND  POINTS  1 AND  2 OF  TRIANGLE 
IF  ITRANSFII  , 1 )- TRANSF  1 1 PT1 , 1 ) ) 360,350,360 
AM  IP  »99999. 

GO  TO  370 

AM 1P*( TRANSFII, 3)- TR AN SFIIPTl*3))/!TRANSF(I»l)- 
TRANSF ! I PT1 , 1 ) ) 

IF  (ABS!AM1P).LT..0001 ) AMiP-0.0 
IF  !TRANSF!I,I)- TRANSF 1 I PT2 , ) ' * 390,380,390 
AM2P  *99999. 

GO  TO  400 

AM2P*I TRANSF < 1 , 3 J-TRANSF 1 1 PT2 ,3 ) ) / < TRANSF ( 1,1)- 
TRANSF1 IPT2.1 ) » 

IF  (ABS(AM2P).LT..000i)  AM2P-0.0 

IAM1P.GE.AM1. AND. AMIP.LE.AM3 I ISL0SW=ISL0SW+1 
(AM2P.LE.AM1.AND.AM2P.GE.AM2)  I SLOSW=  ISLOSW-H 
l ISL0SW.EQ.2)  GO  TO  440 

! TRANSF! 1,1) -TRANSF II PT3,1) ) 420,410,420 
AM  3P  *99999. 

GO  TO  430 

AM3P* ITRANSFII, 3) -TRANSF! IPT3,3) I/! TRANSF!  1,11- 
TRANSF ( I PT3, 1 ) ) 

IF  (ABS(AM3P).LT.. 0001 ) AM3P-0.0 
IF  IAM3P.LE.AM2. AND. AM3P.GE.AM3)  ISLOSW*  ISLOSW-H 
IF  ! I SLOSH. NE. 2 ) GO  TO  460 

IF  1 TRANSF ( I, 2). GE. TRANSF ! I PT1, 2). AND. TRANS  FI  1,2) 

. GE. TRAN SF 1 1 PT2 , 2). AND. TRANSF! 1,21. GE. TRANSF!  IPT3,  2 ) ) 

GO  TO  450 

EQUATION  OF  PLANE  CONTAINING  TRIANGLE 
A*TRANSF ( I PT1 ,2 ) * ( TRANSF 1 1 PT2 , 3) -TRANSF! IPT3.3))- 
TRANSF ( I PT1 , 3 ) * 1 TRANSF ! I PT2 ,2 ) — TRANSF ( I PT3, 2 ) ) ♦ 

( TRANSF 1 1 PT2 ,2 ) *TRANSF ( I PT3 ,3 ) )-! TRANSF ( I PT3»2 )* 

TRANSF ( IPT2 ,3) ) 

B«-ITRANSFIIPTl,l )*( TRANSF ( I PT2, 3) -TRANSF! I PT3,3 ) )- 

TRANSF 1 1 PT1 ,3)* I TRANSF 1 1 PT2 , 1 ) -TRANSF (IPT3,1))4 

TRANSF ! I PT2 , 1 )*TRANSF I IPT3 ,3) -TRANSF! I PT3, 1 ) *T RANSF! IPT2, 

3)1 
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1 

2 

3 

1 

2 

3 

4 


450 

INTEL ( I ) ■! 

460 

CONTINUE 

IF  (IR.EO.O) 

c***» 

CHECK  POINTS  ON  REF 

LL  =NONODE*  LL 

HH=MM*NONODE 

470 

CONTINUE 

480 

CONTINUE 

C "TRANSF  (IPTl»l)*(TRANSF(IPT2,2) -TRANS FI  I PT 3,211- 
TR ANSF ( I PT1 ,2 ) * ( TRANSF ( I PT2  » i ) -TRANSF ( I PT3 , l )) ♦ 

TRANSF ( I PT2, 1 >*TRANSF  I IPT3  ,2 ) -TRANSF ( I PT3 , 1 ) *TRANSF(  IPT2,. 
2) 

0—  ( TRANSF  ( I PT1 , 1 ) *( TRANSF  ( 1 PT2  *2 ) *TRANSF(  I PT3, 3 »- 
TRANSF ( IPT3*2I* TRANSF ( I PT2  *3>  >— TRANSF! I PT1 ,2  >* 

( TRANSF ( I PT2 , 1 ) *TRANSF ( I PT3 ,3 ) -TRANSF ( I PT3, H * 

TRANSF ( 1 PT2  1 3) ) + TRANSF (IPT1»3I*(TRANSF( IPT2.1I* 

TRANSF ( I PT3t 2 l-TRANSF ( I PT3 • 1 ) *TRANSF( t PT2  * 2 1 ) t 
YPLANE*- < A* TRANSF ( I f 1 )+C*TRANSF ( I , 3 )+0) /B 
WRITE  (6,1201  A, B,C,D, YPLANE, TRANSF ( I ,2 », INTEL (II 
IF  (YPLANE.6E. TRANSF (I ,2)1  GO  TO  460 


GO  TO  480 

LECTIVE  TRIANGLES  ALSO. 


RETURN 

END 
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SUBROUTINE  PL0T3D 

C****  THIS  ROUTINE  USES  A MASKING  TECHNIQUE  WHEREIN  THE  UPPER  AND  LOWER  Z'S 
C**«*  FOR  A GIVEN  X ON  THE  PROJECTILE  ARE  STOREO.  LATER,  WHEN  THE  TARGET  IS 
C**»*  PLOTTED,  LINES  OR  PORTIONS  THEREOF  ARE  PLOTTED 
C**+*  WHICH  ARE  ABOVE  OR  BELOW  THIS  MASK  BUT  NOT  WITHIN  IT. 

INCLUDE  COMDIM 
INTEGER  HIGH 
DATA  IFIRST/O/ 


♦♦♦♦♦FORMAT  STATEMENTS •♦♦♦♦ 

IIO  FORMAT! • LOW  * ',16, 'IN  STKATEMENT  I ADD  * lOOOO  - ICHECK, 
♦CHANGE  1CC00  TO  0 IF  LOW  > 20000,  OR  20000  IF  LOW  < 2') 
120  FORMAT  (•  ERR-LOW  MASK  > HIGH  MASK', 3161 


DO  670  K*1,NPNTS 
IF (K  .EQ.  1IL0CSW  * 0 
IX*(XP!K)MOOO.) 

IY*(ZP(K)*1000.) 

IF  (K.NE.l.AND.IMAT.EQ.l)  GO  TO  280 
(LOC)  IS  THE  POSITION  OF  THE  PREVIOUS  POINT  WITH  RESPECT 
TO  THE  MASK 
♦1  ABOVE  THE  MASK 
0 WITHIN  THE  LIMITS  OF  THE  MASK 
-1  BELOW  THE  MASK 

PROCEDURE  FOR  INITIAL  POINT  OF  EACH  LINE 
LOCATE  INITIAL  POINT  WITH  RESPECT  TO  THE  MASK  THEN 
UPDATE  THE  MASK 

IMSKX*IX*SIN (ALPHA ) *1 Y*COS (ALPHA  I 
IF ( IFIRST  .EQ.  1 ) GO  TO  125 
ICHECK  * IMSKX  ♦ IMSKX 
I ADD  * 10C00  - ICHECK 
WRITE (6, 888 1 ICHECK, I ADD 
888  FORM AT( 215) 

IFIRST  * 1 

125  LOW* IMSKX  ♦ IMSKX  ♦ I ADD 

IF  ((LOW. GT.20000.0R.L0W.LT. 2). AND. IMAT.EQ.l)  WRITE  (6,110)  LOW 
C WRITE(6, 178) IX, IY, IMSKX, LOW 

HIGH*L0W-1 

IF  ((HIGH. GE.ISTART. AND. HIGH. LE.IEND). OR. IMAT.EQ.l)  GO  TO  130 
MLOW *-10000 
MH IGH  =— 1 COOO 

I F ( H I GH  .LT.  I START ) LOC SW  * -1 
IF( HIGH  .GT.  IEND ILOCSW  * l 
IF(K  .EQ.  1IL0C0SW  * LOCSW 

I F ( K .NE.  1 .AND.  LOCOLD  .NE.  0 .AND.  LOCSW  .EQ.  LOCOSW) 

♦GO  TO  150 

GO  TO  140 

130  ML  OW*MA  SK ( LOW) 

MHIGH*MASK(HIGH) 

140  CONTINUE 

IF  ( K .NE. 1 .AND. I MAT.EQ.2 ) GO  TO  280 
IF  ( MHIGH-I Y)  180,180,160 
150  CALL  1PL0T  (IX,1Y,2) 

L OC  OL  D *♦  1 
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* 


> I 


non  oo  o oooo 


160 

170 


GO  TO  220 

IF  (MLOW-IY)  170f  200,200 
LOCOLD-O 
GO  TO  210 
180  CONTINUE 

IF  (IMAT.EQ.l)  MA  SK ( H I GH ) ■ I Y 

IF  (ML0W.EQ. -10000. AND. IMAT.EQ.l)  MASK(LOW»«IY 
L OCOL  D ■+ 1 
GO  TO  210 

190  WRITE  (6,120)  MHIGH, JY.MLOW 

STOP 

200  CONTINUE 

IF  (IMAT.EQ.l)  MA  SK (LOW) -I  Y 
LOC  OLD  — 1 

C MOVE  THE  RAISED  PEN  TO  THIS  INITIAL  POINT 

210  IF  (LOCOLD.NE.O. OR. IMAT.EQ.l)  CALL  I PLOT  (IX,IY,3> 

220  JX-IX 

J Y*  I Y 
KX-IX 
KY-I  Y 
GO  TO  670 

SPECIAL  CASE  WHERE  CHANGE  IN  X COORDINATE  IS  ZERO 
A SPECIAL  PROVISION  IS  MADE  AT  THIS  POINT  SO  THAT  A LINE 
WILL  NOT  MASK  ITSELF  AS  LONG  AS  THE  X COORDINATE  REMAINS 
CONSTANT 

280  IF  (IX.NE.JX)  GOTO  290 

JY  = I Y 
GO  TO  310 

COMPUTE  CONSTANTS  FOR  LINEAR  INTERPOLATION 
290  YINC*FLOAT(IY-JY)/ABS(FLOAT(IX-JX)) 

INCX*( IX-JX)/IABS(IX-JX) 

YJ  -J  Y 

PERFORM  LINEAR  INTERPOLATION  AT  EACH  INCREMENTAL  STEP  ON 
THE  X AXIS 
300  JX-JX+INCX 

YJ  *YJ+Y INC 
JY-YJ+.5 

LOCATE  THE  CURRENT  POINT  WITH  RESPECT  TO  THE  MASK  AT  THAT 
POINT  THEN  PLOT  THE  INCREMENT  AS  A FUNCTION  OF  THE 
LOCATION  OF  THE  PREVIOUS  POINT  WITH  RESPECT  TO  ITS  MASK 
310  CONTINUE 

IMSKX = JX*S IN ( ALPHA )*JY*COS( ALPHA) 

LOW-IMSKX  ♦ IMSKX  ♦ IADO 

IF  ( ( L OW.GT. 20000. OR. LOW. LT. 2). AND. I MAT . EQ. 1 ) WRITE  (6.110)  LOW 
HIGH-LOW-1 

IF  ( (HIGH. GE.ISTART. AND. HIGH. LE.IEND). OR. IMAT.EQ.l)  GO  TO  320 
MLOW— 10000 
MHIGH  — 10000 

IF ( HIGH  .LT.  I START) LOC SW  - -1 
IFIHIGH  .GT.  I END ) LOC  SW  - 1 
GO  TO  330 

320  MLOW-MASK ( LOW) 

MHIGH-MASK(HIGH) 

330  IF  (MHIGH- JY)  360.350,340 

340  IF  (MLOW-JY)  370,390,390 
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IF  (MLOW-JY)  360,380,190 
THE  CURRENT  POINT  IS  ABOVE  THE  MASK 
L OC  *♦  1 

IF ( I MAT  .EQ.  1 ) MASK( HIGH ) « JY 

I F ( IMAT  .EQ.  1 .AND.  MLOW  .EQ.  -10000) MASK ( LOW ) * JY 
I F ( IMAT  .EQ.  1 )G0  TO  650 

IF  ( LOCSW. NE.O. AND. LOCO LD.NE.O)  GO  TO  640 
IF  (LOCOLO)  470,450,590 
THE  CURRENT  POINT  IS  WITHIN  THE  MASK 


370  L OC  =0 

LOCSW  * 0 ) 

I F ( IMAT  .EQ.  1 )G0  TO  650 
IF  (LOCOLD)  470,650,470 

380  IF  (LOCOLO)  390,370,360 

C THE  CURRENT  POINT  IS  BELOW  THE  MASK 

390  L OC*- 1 

IF( IMAT  .EQ.  l)MASK(LOW)  * JY 
IFl IMAT  .EQ.  1 ) GO  TO  650 

IF  (LOCSW. NE.O. AND. LOCOLD. NE.O)  GO  TO  640 
IF  (LOCOLD)  640,460,470 
400  IF  (LOCOLD)  420,650,410 

C****  pL0T  FRqm  ABOVE  MASK  TO  TOP  EDGE  OF  MASK 
410  I ARG1 =HIGHY 

IARG2*MHIGH 
GO  TO  660 

C****  plot  FROM  BELOW  MASK  TO  BOTTOM  EDGE  OF  MASK 
420  IARG1=L0WY 

I ARG2=ML0W 
GO  TO  660 

430  IF  (ABS(KY-IARG1).GT.ABS( JY-IARG2))  GO  TO  440 

IF  ( IARG1.NE. -10000)  CALL  IPLOT  (KX,IARG1,2) 

GO  TO  670 

440  IF  ( IARG2.NE. -10000)  CALL  IPLOT  (JX,IARG2,2I 

GO  TO  670 

450  CALL  IPLOT  (JX,JY,3) 

GO  TO  590 

460  CALL  IPLOT  (JX,JY,3) 

GO  TO  640 

470  IF  (LOCOLD. NE.O. OR. IMAT. EQ.l)  CALL  IPLOT  (KX,KY,2) 

IMSKX =KX* SIN ( ALPHA )*KY*COS (ALPHA) 

MSKL  * I MSKX  ♦ IMSKX  ♦ I ADD 

IF  ( (MSKL. GT.20000.0R.MSKL.LT. 2) .AND. IMAT. EQ.l)  WRITE  (6,110) 

1 MSKL 

MSKH*MSKL-1 

IF  ((MSKH.GE.ISTART.ANO.MSKH.LE.IEND). OR. IMAT. EQ.l)  GO  TO  480 
LOWY *-10000 
HIGHY *-10000 

IFtMSKH  .LT.  I START) LOCSW  * -1 
IF(MSKH  .GT.  IEND ) LOCSW  * 1 

GO  TO  490  , 

480  LOWY -MASK ( MSKL ) 

HIGHY -MASK(MSKH) 

490  IF  (LOC)  500,400,510 

500  CONTINUE 

GO  TO  520 
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510  CONTINUE 

IF  ( IMAT.EQ. 1 ) MA SK ( MSKL ) »KY 
520  IF  (LOCOLD)  540,530,550 

530  IF  (LOC)  610,650,600 

C PLOT  FROM  BELOW  MASK  TO  ABOVE  MASK 

540  IF  (LOCSW.NE.O)  GO  TO  640 

I ARG1 =LOWY 
IARG2*MHIGH 
IARG3=HIGHY 
IARG4=ML0W 
ASSIGN  590  TO  N 
GO  TO  560 

C PLOT  FROM  ABOVE  MASK  TO  BELOW  MASK 

550  IF  (LOCSW.NE.O)  GO  TO  640 

I ARG1 =HIGHY 
I ARG2=ML0W 
I ARG3=L0WY 
IARG4=MHIGH 
ASSIGN  640  TO  N 

560  IF  ( ABS(KY-IARGl) .GT.ABS( JY-IARG2) ) GO  TO  580 

IF  ( IARG1.NE. -10000)  CALL  IPLOT  <KX,IARG1,2) 

C GET  TO  OTHER  SIDE  OF  MASK 

IF  ( IARG3.NE. -10000)  CALL  IPLOT  (KX.IARG3.3) 

570  CALL  IPLOT  (JX,JY,2) 

GO  TO  N 

580  IF  ( IARG4.NE. -10000)  CALL  IPLOT  (JX,IARG4,2) 

C**»*  GET  TO  OTHER  SIDE  OF  MASK 

IF  ( IARG2.NE. -10000)  CALL  IPLOT  (JX,IARG2,3> 

GO  TO  570 

C PLOT  FROM  ABOVE  THE  MASK  TO  ABOVE  THE  MASK 

590  CONTINUE 

GO  TO  650 

C PLOT  FROM  INSIDE  MASK  TO  ABOVE  MASK 

600  I ARG1=HIGHY 

I ARG2=MHIGH 
ASSIGN  590  TO  N 
GO  TO  620 

C PLOT  FROM  INSIDE  MASK  TO  BELOW  MASK 

610  IARG1=L0WY 

I ARG2=ML0W 
ASSIGN  640  TO  N 

620  IF  (ABS(KY-IARGl).GT.ABS(JY-I ARG2 ) ) GO  TO  630 

IF  ( IARG1.NE. -10000)  CALL  IPLOT  (KX,IARG1,3) 

GO  TO  570 

630  IF  ( IARG2.NE. -10000)  CALL  IPLOT  (JX,IARG2,3) 

GO  TO  570 

C PLOT  FROM  BELOW  THE  MASK  TO  BELOW  THE  MASK 

640  CONTINUE 

L 0CSW*0 
650  KX*JX 

KY«JY 

660  L0C0LD«L0C 

IF  (JX.NE.IX)  GO  TO  300 

IF  (LOC. NE.O. OR. IMAT.EQ. 1)  CALL  IPLOT  (JX,JY,2) 
IF  (LOCOLD. EQ.O. AND. LOC. EQ.O)  GO  TO  670 
IF  (LOC. EQ.O)  GO  TO  430 
670  CONTINUE 

RETURN 
END 
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SUBROUTINE  I PLOT! I X, I Y, I P) 
INCLUDE  COMDIH 
2 AX  * IX/IOOO. 

AY  * IY/ICOO. 

AX  * (AX  - XP(5II/XP<6) 

AY  * (AY  - ZP(5>)/ZP(6) 
CALL  PLOT( AX t AY 1 1 P ) 

RETURN 

END 


) 


J 

i 
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r 
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SUBROUTINE  CKMASK 

C***»  THIS  ROUTINE  CHECKS  MASK  LIMITS  TO  SEE  IF  THEY  FIT  WITHIN  RANGE  1-20001 
INCLUDE  COMDIM.LIST 


♦♦•♦♦FORMAT  STATEMENTS 

110  FORMAT  (•  MASK  COMPLETELY  EMPTY* I 
120  FORMAT  (2I6.F10.3) 

C 

DO  130  I >1 1 20000*2 
IF  (MASK! I ).NE. -10000)  GO  TO  140 
130  CONTINUE 

WRITE  (6.110) 

ISTART-30000 
I END-30000 
GO  TO  170 
140  ISTART-I 

DO  150  I-ISTART .20000.2 
IF  (MASK(I).EQ. -10000)  GO  TO  160 
150  CONTINUE 

160  I END-I-2 

WRITE  (6,120)  ISTART.IEND, ALPHA 
170  CONTINUE 
RETURN 
END 


APPENDIX  E 


RETRIEVAL  RUNSTREAM 


IASG,A  RDUM2S • 

•use  22,RDunaa. 

•MAP  .TPFf.ABS 
IN  RLIB.RETREU 
LIB  RUB* 

LIB  SEAPPD#TEKLIB2. 
IXQT 
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APPENDIX  F 


PROCEDURE  TO  DETERMINE  IF  A POINT  IS  HIDDEN  BY  A TRIANGLE 


The  following  procedure  is  used  in  the  program  to  determine  whether 
a point  (P)  is  hidden  by  a triangle  defined  by  points  A,  B,  C. 

1.  P is  not  hidden  by  AABC  if  one  of  the  following  conditions  is 

met: 


2.  Arrange  the  points  defining  the  triangle  in  counterclockwise 
order  (PI,  P2,  P3)  such  that  Zj  = min  (zA,  zB,  zc) . (If  zm  = zn  = Zj 
then  PI  is  point  wherein  = min  (xm,  x^).  Let  the  slope  between  any 
two  points  A and  B be  M(A,  B) . Then  if  2 of  the  following  3 conditions 
are  met  the  point  (P)  lies  within  the  area  of  the  triangle  and  must  be 
further  checked  (see  Step  3).  Otherwise  it  is  not  hidden  by  AABC. 


Determine  the  equation  of  the  plane  defined  by  PI,  P2,  P3. 
(at  xp,  zD)  < yp  then  the  point  (P)  is  hidden  by  the  triangle 


DISTRIBUTION  LIST 


No.  of  No.  of 

Copies  Organization  Copies  Organization 


12  Commander 

Defense  Documentation  Center 
ATTN:  DDC-TCA 
Cameron  Station 
Alexandria,  VA  22314 


1 


Commander 

US  Army  Missile  Materiel 
Readiness  Command 
ATTN:  DRSMI-AOM 
Redstone  Arsenal,  AL  35809 


1 


Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCDMD-ST,  N.  Klein 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 


1 Commander 

US  Army  Aviation  Research 
and  Development  Command 
ATTN:  DRSAV-E 
P.  0.  Box  209 
St.  Louis,  M0  63166 


1 


Director 

US  Army  Air  Mobility  Research 
and  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 


1 


Commander 

US  Army  Electronics  Research 
and  Development  Command 
Technical  Support  Activity 
ATTN:  DELSD-L 
Fort  Monmouth,  NJ  07703 


1 


Commander 

US  Army  Communications  Rsch 
and  Development  Command 
ATTN:  DRDCO-SGS 
Fort  Monmouth,  NJ  07703 


1 Commander 

US  Army  Missile  Research 
and  Development  Command 
ATTN:  DRDMI-R 
Redstone  Arsenal , AL  35809 


1 Commander 

US  Army  Tank  Automotive 
Research  § Development  Cmd 
ATTN : DRDTA-UL 
Warren,  MI  48090 

2 Commander 

US  Army  Armament  Research 
and  Development  Command 
ATTN:  DRDAR-TSS  (2  cys) 

Dover,  NJ  07801 

1 Commander 

US  Army  Armament  Materiel 
Readiness  Command 
ATTN:  DRSAR-LEP-L,  Tech  Lib 
Rock  Island,  IL  61299 

5 Commander 

US  Army  Materials  and 

Mechanics  Research  Center 
ATTN:  DRXMR-T,  Mr.  J.  Bluhm 

DRXMR-T,  Dr.  D.  Roylance 
DRXMR-T,  Dr.  A. F.  Wilde 
Dr.  J.  Mescal 1 
DRXMR-ATL 

Watertown,  MA  02172 
1 Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SL,  Tech  Lib 
White  Sands  Missile  Range 
NM  88002 


79 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 

1 Commander 

US  Army  BMD  Advanced 
Technology  Center 
ATTN:  Mr.  P.  Boyd,  BMDATC-M 
P.  0.  Box  1500 
Huntsville,  AL  35807 

3 Commander 

Naval  Weapons  Center 
ATTN:  Code  4057 

Code  4011,  E.Lundstrom 
Code  3813,  M.  Backman 
China  Lake,  CA  93555 

2 ADTC/DLJW  (MAJ  D.  Matuska, 

LTC  J.  Osborn) 

Eglin  AFB,  FL  32542 

3 Honeywell,  Inc. 

Government  § Aerospace 

Products  Division 
ATTN:  Mr.  J.  Blackburn 
Dr.  G.  Johnson 
Mr.  R.  Simpson 
600  Second  Street,  NE 
Hopkins,  MN  55343 


No.  of 

Copies  Organization 

1 Kaman  Sciences  Corporation 
ATTN : Dr . P . Snow 
1500  Garden  of  the  Gods  Road 
Colorado  Springs,  CO  80933 

3 Sandia  Laboratories 
ATTN : Dr . W . Herrmann 
Dr.  L.  Bertholf 
Dr.  J.  W.  Nunziato 
Albuquerque,  NM  87115 

1 Drexel  University 

Dept  of  Mechanical  Engineering 
ATTN:  Dr.  P.  C.  Chou 
32nd  and  Chestnut  Streets 
Philadelphia,  PA  19104 

T Stanford  Research  Institute 
ATTN:  Dr.  L.  Seaman 
Dr . D . Curran 
Dr.  D.  Shockey 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 


Aberdeen  Proving  Ground 
Dir,  USAMSAA 
Cdr,  USATECOM 

ATTN:  DRSTE-SG-H 


